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ABSTRACT 



I Aims. We have performed a comprehensive multiwavelength analysis of a sample of 20 starburst galaxies that show a 

substantial population of very young massive stars, most of them classified as Wolf-Rayet (WR) galaxies. In this paper, 
the forth of the series, we present the global analysis of the derived photometric and chemical properties. 
Methods. We compare optical/7V/7? colours and the physical properties (reddening coefficient, equivalent widths of 
the emission and underlying absorption lines, ionization degree, electron density, and electron temperature) and 
chemical properties (oxygen abundances and N/0, S/0, Ne/0, Ar/0, and Fe/0 ratios) with previous observations 
and galaxy evolution models. We compile 41 independent star-forming regions -with oxygen abundances between 
12-|-log(0/II)= 7.58 and 8.75-, of which 31 have a direct estimate of the electron temperature of the ionized gas. 
Results. According to their absolute B-magnitude, many of them are not dwarf galaxies, but they should be during 
, their quiescent phase. We found that both c(II/3) and Wats increase with increasing metallicity. The differences in the 

N/0 ratio is explained assuming differences in the star formation histories. We detected a high N/O ratio in objects 
showing strong WR features (HCG 31 AC, UM 420, IRAS 0828-1-2816, III Zw 107, ESO 566-8 and NGC 5253). The 
ejecta of the WR stars may be the origin of the N enrichment in these galaxies. We compared the abundances provided 
^ , by the direct method with those obtained through empirical calibrations, finding that (i) the Pilyugin method is the 

■ best suited empirical calibration for these star-forming galaxies, (ii) the relations provided by Pettini & Pagel (2004) 
^ ' give acceptable results for objects with 12-|-log(O/H)>8.0, and (iii) the results provided by empirical calibrations based 

I on photoionization models are systematically 0.2 - 0.3 dex higher than the values derived from the direct method. The 

O and N abundances and the N/O ratios are clearly related to the optical/iV/7? luminosity; the dispersion of the data 
is a consequence of the differences in the star-formation histories. The L-Z relations tend to be tighter when using NIR 
04 ' luminosities, which facilitates distinguishing tidal dwarf galaxies candidates and pre-existing dwarf objects. Galaxies 

^ ' with redder colours tend to have higher oxygen and nitrogen abundances. 

, Conclusions. Our detailed analysis is fundamental to understand the nature of galaxies that show strong starbursts, as 

04 ■ well as to know their star formation history and the relationships with the environment. This study is complementary 

' -but usually more powerful- to the less detailed analysis of large galaxy samples that are very common nowadays. 

O ' 

• ' Key words, galaxies: starburst — galaxies: interactions — galaxies: dwarf — galaxies: abundances — galaxies: kinematics 

\l ' and dynamics — stars: Wolf-Rayet 

o 

1. Introduction key ingredient for modelling galaxy properties, because it 

1^ ■ , , , , „ , , . , . . r 1 . . determines f/F, optical and A'^/iZ colours at a given age (i.e., 

■ The knowledge of the chemical composition of galaxies m Lgit^^^rer et al. 1999), nucleosynthetic yields (e.g., Woosley 
^ . particular of dwarf galaxies is vital for understanding their ^ ^^^^^^ ^995^^ dust-to-gas ratio (e.g., Hirashita et al 

^ . evolution, star formation history, stellar nuc eosynthesis, 200I), the shape of the interstellar extinction curve (e.g., 

... the importance of gas mflow and outflow and the enrich- pipya net al 2006), or even the properties of the Wolf-Rayet 

ment ot the mtergalactic medium. Indeed, metaUicity is a g^^rs (|Crowtheij|2007t) 

Send offprint requests to: Angel R. Lopez-Sanchez, e-mail: The most robust method to derive the metallicity in 

Angel.Lopez-Saiichez@csiro.au star-forming and starburst galaxies is via the estimate 

* Based on observations made with NOT (Nordic Optical of metal abundances and abundance ratios, in particu- 

Telescope), INT (Isaac Newton Telescope) and WHT (WiUiam i^r through the determination of the gas-phase oxygen 

Herschel Telescope) operated on the island of La Palma jointly abundance and the nitrogen-to-oxygen ratio. The relation- 

by Denmark, Finland, Iceland, Norway and Sweden (NOT) 1 • 1 , j. j. n- j ^1 1 

or the Isaac Newton Group (INT, WHT) in the Spanish ships between current metalhcity and other galaxy param- 

Observatorio del Roque de Los Muchachos of the Instituto de ^ters, such as colours, luminosity, neutral gas content, star- 

Astrofi'sica de Canarias. Based on observations made at the formation rate, extinction or total mass, constrain galaxy- 

Centro Astronomico Hispano Aleman (CAHA) at Calar Aho, evolution models and give clues about the current stage 

operated by the Max-Planck Institut fiir Astronomie and the of a galaxy. For example, is still debated whether massive 

Instituto de Astrofi'sica de Andaluci'a (CSIC). star formation results in the instantaneous enrichment of 
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the interstellar medium of a dwarf galaxy, or if the bulk 
of the newly synthesized heavy elements must cool before 
becoming part of the interstellar medium (ISM) that even- 
tually will form the next generation of stars. Accurate oxy- 
gen abundance measurements of several H ii regions within 
a dwarf galaxy will increase the understanding of its chem- 
ical enrichment and mixing of enriched material. The anal- 
ysis of the kinematics of the ionized gas will also help to 
understand the dynamic stage of galaxies and reveal recent 
interaction features. Furthermore, detailed analyses of star- 
burst galaxies in the nearby Universe are fundamental to 
interpret the observations of hig h-z star f orming galaxies, 
such as Lyman Break Galaxies ()Erb et al.|[2003f ) , as well as 
quantify the importance of interactions in the triggering of 
the star-formation bursts, which seem to be very common 
at higher redshifts (i.e., Kauffmann & White 1993; Springer 
et al. 2005). 

The comparison of the metallicity (which reflects the gas 
reprocessed by stars and any exchange of gas between the 
galaxy and its environment) with the stellar mass (which 
reflects the amount of gas locked up into stars) provides 
key clues about galaxy formation and evolution. These 
analyses have shown a clear correlation between mass and 
metallicity. In practice, luminosity has been used as sub- 
stitute of mass because of the difficulty of deriving reli- 
able galaxy masses, yielding to the so-called metallicity- 
luminosity relation (i.e., Rubin et al. 1984; Richer & McCall 
1995; Salzer et al. 2005), although in recent years mass- 
metallicity relations are also explored (i.e., Tremonti et al. 
2004; Kewley & Elisson, 2008), and are studied even at high 
redshifts (i.e., Kobulnicky et al. 1999; Pettini et al. 2001; 
Kobulnicky & Kewley 2004; Erb et al. 2006; Liang et al. 
2006) . The evolution of such relationships are now predicted 
by semi-analytic models of galaxy formation within the A- 
cold dark matter framework that include chemical hydrody- 
namic simulations (De Lucia et al. 2004; Tissera et al. 2005; 
De Rossi et al. 2006; Dave & Oppenheimer 2007). Ironically, 
today the main problem is not to estimate the mass of a 
galaxy but its real metallicity, so that different methods in- 
volving direct estimates of the oxygen abundance, empirical 
calibrations using bright emission-line ratios or theoretical 
methods based on photoionization models yield very differ- 
ent values (i.e., Yin et al. 2007; Kewley & Elisson, 2008). 

Hence precise photometric and spectroscopic data, in- 
cluding a detailed analysis of each particular galaxy that 
allows conclusions about its nature, are crucial to address 
these issues. We performed such a detailed photometric 
and spectroscopic study in a sample of strong star-forming 
galaxies, many of them previously classified as dwarf galax- 
ies. The majority of these objects are Wolf-Rayet (WR) 
galaxies, a very inhomogeneous class of star-forming ob- 
jects which share at least an ongoing or recent star for- 
mation event that has produc ed stars sufficiently massive 
to ev olve into the WR stage (jSchaerer. Contini fc Pindaol 
Il999|) . However, WR features in the spectra of a galaxy 
provides useful information about the star-formation pro- 
cesses in the system. As the first WR stars typically ap- 
pear around 2-3 Myr after the starburst is initiated an d 
disappear within some 5 Myr ()Mevnet &: MaedeH l2005l) , 
their detection gives indications about both the youth and 
strength of the burst, offering the opportunity to study 
an approximately coeva l sample of very young starbursts 
(jSchaerer fc Vaccal[l998l ). 



The main aim of our study of the formation of mas- 
sive stars in starburst galaxies and the role that the in- 
teractions with or between dwarf galaxies and/or low sur- 
face bright ness objects have in its triggerin g mechanism. 
In Paper I (iLopez-Sanchez fc EstebanI I2008D we described 
the motivation of this work, compiled the list of the 20 
analysed WR galaxies (Table 1 of Paper I), the majority 
of them showing several sub-regions or objects within or 
surrounding them, and presented the results of the op- 
tical/ NIR broad-band and Ha photometry. In Paper II 
(jLopez- Sanchez fc EstebanI I2009D we presented the results 
of the analysis of the intermediate resolution long-slit spec- 
troscopy of 16 WR galaxies of our sample - the results 
for the other four galaxies were published separately. In 
many cases, two or more slit positions were used to analyse 
the most interesting zones, knots or morphological struc- 
tures belonging to each galaxy or even surrou nding ob- 
jects. Paper HI ()L6pez-Sanchez fc Estebanll2010D presented 
the analysis of the O and WR stellar populations within 
these galaxies. In this paper, the forth of the series, we 
globally compile and analyse the optical/ A/i? photomet- 
ric data (Sect. 2) and study the physical (Sect. 3) and 
chemical (Sect. 4) properties of the ionized gas within our 
galaxy sample. Thirty-one up to 41 regions have a direct 
estimate of the electron temperature of the ionized gas, and 
hence the element abundances were derived with the direct 
method. Section 4 includes the analysis of the N/0 ratio 
with the oxygen abundance, a discussion of the nitrogen 
enrichment in WR galaxies, a study of the a-elements to 
oxygen ratio with the oxygen abundance, and the compar- 
ison of the results provided by the most common empirical 
calibrations with those derived following the direct method 
(the Appendix compiles all metallicity calibrations used in 
this work) . Section 5 analyses the metallicity-luminosity re- 
lations obtained with our data. Section 6 discusses the rela- 
tions between the metallicity and the optical/ NIR colours. 
Finally, we list our main conclusions in Sect. 7. 

The final paper of the series (Paper V) will compile 
the properties derived with data from other wavelengths 
(UV, FIR, radio, and X-ray) and complete a global analysis 
of all available multiwavelength data of our WR galaxy 
sample. We have produced the most comprehensive data set 
of these galaxies so far, involving multiwavelength results 
and analysed according to the same procedures. 



2. Global analysis of magnitudes and colours 

Our optical and NIR broad-band photometric results for 
the galaxy sample were presented in Paper 1. These data 
allowed us the analysis of the optical and NIR magnitudes 
and the colours of the galaxies and surrounding dwarf ob- 
jects. Table [T] compiles the optical/A'/i? results for the in- 
dividual galaxies, not considering regions within them or 
nearby diffuse objects. This table shows the colour excess, 
E{B—V) (derived with the Balmer decrement in our optical 
spectra, see Paper II), the absolute S-magnitude (both cor- 
rected, Mb, and uncorrected, M^'^ , for extinction), all the 
optical/A/i? colours, and the age of the most recent star- 
formation burst (the young population, derived from our 
Ha images) and the minimum age of the old stellar popula- 
tion (usually estimated from the low luminosity component 
or regions without nebular emission using our optical/A/i? 
broad-band images). 
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Table 1. Compilation of the broad-band photometric data" for the individual galaxies analysed in this work. 



Galaxy 


E{B-V) M^^ 


Mb 


U -B 


B-V 


V -R 


V -J 


J-H 


H-Ks 


B. Age 


UC Age 




(a) 


(b) 




id) 




(a) 


id) 










id) 






HCG 31 AC 


0.06 


-19.18 


-19.43 


-0.60+0.06 


0.03+0.08 


0.12+0.08 


0.20+0.10 


0.13+0.10 


0, 


.15+0.12 


5.0 


100 


HCG 31 B 


0.18 


-17.96 


-18.71 


-0.38+0.08 


0.17+0.06 


0.06+0.06 


0.14+0.10 


0.13+0.10 


0, 


.12+0.10 


7.0 


100 


HCG 31 E 


0.06 


-15.51 


-15.76 


-0.65+0.10 


-0.03+0.10 


0.20+0.09 


0.29+0.12 


0.05+0.10 


0, 


.18+0.12 


6.0 


- 


HCG 31 Fl 


0.20 


-14.93 


-15.76 


-0.99+0.12 


-0.07+0.12 


-0.04+0.10 


-0.17+0.14 


0.04+0.17 


0, 


.29+0.30 


2.5 





HCG 31 F2 


0.09 


-13.97 


-14.34 


-1.01+0.12 


-0.09+0.12 


-0.02+0.10 


0.01+0.16 


0.08+0.30 


0, 


.20+0.50 


2.5 





HCG 31 G 


0.06 


-18.63 


-18.88 


-0.43+0.09 


-0.01+0.08 


n 1 /I -Un 


.Uo 


0.45+0.08 


0.12+0.10 


0, 


.13+0.10 


6.0 


100 


Mkn 1087 


0.17 


-21.45 


-22.15 


-0.41+0.08 


0.17+0.08 


n ori-un 
U.zUitU. 


.Uo 


0.52+0.06 


0.20+0.06 


0, 


.16+0.06 


6.0 


100 


Mkn 1087 N 


0.10 


-17.65 


-18.06 






-0.05+0.06 


n 1 /I _i_ri 


1 n 

.lU 


0.21+0.08 


0.18+0.08 


0, 


.13+0.08 


7.0 


— 


Mkn 1087 #1 


0.079 


-16.05 


-16.34 


-0.75+0.15 


-0.01+0.10 


n 1 ri-un 

U.lUltU. 


.Uo 










6.0 


- 


Mkn 1087 #3 


0.07" 


-16.91 


-17.20 


0.08+0.30 


0.11+0.06 


n op 1 n 


.Uo 


0.64+0.10 


0.50+0.20 








150 


Haro 15 


0.11 


-20.41 


-20.87 


-0.52+0.08 


0.26+0.08 


n Qo_i_n 


.Uo 


0.17+0.08 


0.58+0.08 


0, 


.22+0.08 


5.0 


500 


Mkn 1199 


0.15 


-20.06 


-20.68 


-0.44+0.06 


0.46+0.06 


U.zyitU. 


.UD 


1.30+0.07 


0.55+0.08 


0, 


.34+0.08 


8.0 


500 


Mkn 1199 NE 


0.11 


-17.11 


-17.57 


0.16+0.08 


0.51+0.08 


n Q/i_i_n 


.Uo 


1.29+0.08 


0.62+0.10 


0, 


.20+0.10 


12.0 


500 


Mkn 5 


0.20 


-14.74 


-15.57 


-0.41+0, 


,06 


0.44+0.06 




.UD 


0.81: 


0.52+0.03 


0, 


.38+0.04 


5.0 


500 


IRAS 08208-1-2816 


0.17 


-20.59 


-21.29 


-0.49+0, 


,06 


0.22+0.06 




.Uo 


1.03+0.08 


0.54+0.08 


0, 


.22+0.10 


5.5 


500 


IRAS 08339-1-6517 


0.16 


-20.91 


-21.57 


-0.51+0, 


.08 


0.01+0.08 


U.zDitU. 


.Uo 


1.36+0.06 


0.64+0.05 


0, 


.23+0.06 


4.5 


150 


IRAS 08339-1-6517 C 


0.13 


-17.67 


-18.21 


-0.16+0, 


,10 


0.20+0.08 


U.zDltU. 


.Uo 


1.45+0.12 


0.21+0.25 


0, 


.68+0.28 


5.5 


250 


POX 4 


0.06 


-18.54 


-18.79 


-0.68+0, 


,03 


0.29+0.02 


n '^9-1-0 


.U^ 


0.42+0.08 


0.28+0.08 


0, 


.15+0.10 


3.5 


250 


POX 4 Comp 


0.12 


-14.86 


-15.36 


-0.02+0, 


.06 


0.25+0.02 


0.30+0, 


.04 


0.87+0.10 


0.7: 




0.3: 


5.0 


300 


UM 420 


0.06 


-19.30 


-19.55 


-0.80+0, 


,06 


0.31+0.06 


0.13+0, 


,06 


0.77+0.12 


0.41+0.12 


0, 


.12+0.16 


4.5 


200 


BBS 0926+606 A 


0.08 


-16.96 


-17.29 


-0.75+0, 


,06 


0.01+0.06 


0.14+0, 


,06 


0.54+0.06 


0.21+0.06 


0, 


.15+0.08 


4.8 


200 


BBS 09264-606 B 


0.12 


-16.87 


-17.37 


-0.51+0, 


.08 


0.08+0.06 


0.20+0, 


.06 


0.83+0.06 


0.29+0.06 


0, 


.18+0.08 


6.7 


100 


BBS 0948-1-532 


0.24 


-17.44 


-18.43 


-1.20+0, 


,06 


-0.12+0.06 


0.16+0, 


,06 










4.6 


100 


BBS 1054-1-365 


0.02 


-13.98 


-14.06 


-0.34+0, 


,06 


0.33+0.06 






0.92+0.08 


0.38+0.12 


0, 


.16+0.15 


4.9 


500 


BBS 1211+540 


0.08 


-12.94 


-13.27 


-0.61+0, 


.06 


0.04+0.06 


0.21+0, 


.06 










4.7 


500 


BBS 1319+579 


0.02 


-18.45 


-18.53 


-0.39+0, 


.06 


0.34+0.06 


0.19+0, 


.06 


1.03+0.08 


0.39+0.12 


0, 


.16+0.20 


3.7 


300 


BBS 1415+437 


0.13 


-14.09 


-14.52 


-0.47+0, 


,06 


0.21+0.06 


0.27+0, 


,06 


0.98+0.08 


0.35+0.10 




0.15: 


3.6 


250 


TTT ^T.T 1 (\'7 
ill ZjW iU( 


0.21 


-19.27 


on 1 A 
-ZD. 14 


-0.42+0, 


.06 


n 1 /I _l_n nft 
U.14itU.U0 


0.22+0, 


.06 


V.oZzLV.lZ 


n A^j-(\ on 
U.4 / ±U.^U 


0, 


.35+0.20 


0.0 


OUU 


Tol 9 


0.31 


-17.98 


-19.26 


-0.34+0, 


.06 


0.24+0.06 


0.22+0, 


.06 


0.83+0.08 


0.68+0.10 


0, 


.27+0.12 


5.8 


500 


Tol 1457-262 Objl 


0.16 


-19.07 


-19.73 


-0.56+0, 


,06 


0.23+0.06 


0.26+0, 


,06 


0.60+0.10 


0.51+0.12 


0, 


.22+0.12 


4.6 


500 


Tol 1457-262 Obj2 


0.169 


-18.31 


-18.97 


-0.42+0, 


.06 


0.34+0.06 


0.36+0, 


.06 


0.90+0.10 


0.58+0.12 


0, 


.27+0.14 


5.2 


500 


Tol 1457-262 #15 


0.169 


-15.82 


-16.48 


-0.25+0, 


.10 


0.39+0.08 


0.39+0, 


.06 


1.10+0.20 








6.4 


400 


Tol 1457-262 #16 


0.169 


-14.03 


-14.69 


-0.10+0, 


,15 


0.45+0.10 


0.40+0, 


,06 










7.0 


500 


ESO 566-8 


0.34 


-19.47 


-20.88 


-0.48+0, 


.06 


0.31+0.06 


0.19+0, 


.06 


1.10+0.10 


0.60+0.12 


0, 


.38+0.14 


4.2 


500 


ESQ 566-7 


0.16 


-18.69 


-19.35 


-0.21+0, 


.08 


0.49+0.06 


0.31+0, 


.06 


1.20+0.10 


0.71+0.16 


0, 


.34+0.16 


4.2 


500 


NGC 5253*^ 


0.17 


-17.23 


-17.92 


-0.41+0, 


.02 


0.27+0.02 


0.21+0, 


.02 


0.81+0.03 


0.53+0.04 


0, 


.19+0.05 


3.5 


300 



° Colour excess, E{B — V), derived from our estimates of the reddening coefficient and assuming iJv'=3.1, EiB — y)=0.692c(H/3). 
^ Absolute B-magnitude, not corrected for extinction. 

" Extinction-corrected absolute B-magnitude, assuming Mb = M^'~^ - Ab = M^'^ - 3.1 X 1.337E(B - V) = M^'^ - 2.868 c(H/3). 

All colours have been corrected for both extinction and emission of the gas, see lLopez-Sanchez fc EstebanI l l2008f) . 
" Age of the most recent star- forming burst (derived using the Ha equivalent width), in Myr. 

f Minimum age of the underlying stellar population (derived via the analysis of the low-luminosity component of the galaxy), in Myr. 
9 EiB — V) was estimated only considering the extinction of the Milky Way. 

^ Optical and NIR magnit udes extracted from the NED and corrected for extinction using an average value of c(H/3)=0.24, see 
ILopez-Banchez et all 1120071 ). 



Our first result from Table [T] is that the actual number 
of dwarf galaxies, defined as Mb > — 18, is not as high as we 
had expected considering the selection criteria of our WR 
galaxy sample. There are two reasons for this: (i) on the one 
hand, the determination of the magnitudes was performed 
in a more accurate way. As our images are deeper than 
those previously obtained, the integrated magnitude of a 
diffuse object is lower than that estimated before, (ii) On 
the other hand, we corrected all our data for extinction, but 
not only considering the effect of the dust in the Milky Way 
as it is usually done, but taking into account the internal 
extinction derived from our spectroscopic data. That is why 
only six galaxies (Mkn 5, SBS 0926+606, SBS 1054+365, 
SBS 1211+540, SBS 1415+437 and NGC 5253) are strictly 
classified as dwarf galaxies following the above defini- 
tion. POX 4, SBS 0948+532 and SBS 1319+579 could 
be also considered dwarf galaxies because Mb > —19. 
Table [1] also lists some tidal dwarf galaxy (TDG) can- 



didates (HCG 31 E, Fl and F2; Mkn 1087 #1 and #3, 
POX 4 Comp) and nearby external objects (Mkn 1087 N, 
Mkn 1199 NE, IRAS 08309+6517 C, Tol 1457-262 #15 and 
#16) surrounding a main galaxy. However, as we remarked 
i n Paper I and in the analysis of the HCG 31 members 
()L6pez- Sanchez. Esteban fc Rodn'gue^ l2004aD . we must 
keep in mind that the i?-magnitude of a starburst is in- 
creased by several magnitudes during the first 10 Myr with 
respect to its brightness in the quiescent phase, so we should 
expect that some of the objects with — 18 > Mb > —19.5 
are indeed defined as dwarf objects during their quiescent 
phase. 

In order to quantify the effect of the correction for ex- 
tinction, we plot the uncorrected absolute B-magnitude 
versus the absolute B-magnitude corrected for extinction in 
Fig. [T] We did not consider the correction for the emission 
of the gas in the absolute magnitude because (i) it is small 
in the i?-filter, less than 0.10 magnitudes and, more impor- 
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Mg (corrected for extinction) 

Fig. 1. Comparison between the reddening-corrected absolute 
_B-magnitudes {Mb) and the uncorrected ones (Mf ^) for all 
individual galaxies analysed in this work. A linear fit to the 
data is shown. 

tant, (ii) we are considering the magnitude of the galaxy as 
a whole, taking into account both the star-forming bursts 
and regions dominated by older stellar populations that 
do not possess any nebular emission. From Fig. [TJ we see 
that the magnitudes corrected for extinction are on aver- 
age around 0.60 magnitudes lower than when this effect is 
not considered. As all data lie in a narrow band, we per- 
formed a linear fit, finding the following relation between 
both magnitudes: 

Ms = 1.02M]^'=-0.19. (1) 

For Mb = —18, the magnitude difference is AMb ~0.56. 
Transforming this value to luminosity, the consideration of 
the correction for extinction means that one has to multiply 
the observed B-luminosity of a galaxy by a factor between 

1.6 (for Mb 16) and 1.8 (for Mb ~ -22). We note 

that there is a slight dependence on the extinction with the 
absolute magnitude of the galaxy, that is, the correction 
for extinction is higher at lower absolute magnitudes. This 
suggests a higher absorption of the light in brighter systems 
(more amount of dust). We will get the same result when 
we analyse the relation between the reddening coefficient 
and the warm dust mass (Paper V). 

As we explained in Paper I, we compared our 
optical/NIR colours (corrected for extinction and 
emission of the ionized gas) with the predictions 
given by three different population synthesis models, 
STARBURST99 jLc ithcrcr ct_al. 1 999j), PEGASE.2 
( Fioc &: Rocca-Volmerangel 119971 ) and iBruzual fc Charloti 
(2003), to estimate the age of the dominant stellar pop- 
ulation of the galaxies, the star-forming regions, and the 
underlying stellar component. We assumed an instanta- 
neous burst with a Salpeter IMF a total mass of 10® Mq, 
and a metallicity of Z/Zq= 0.2, 0.4 and 1 (chosen as a 
function of the oxygen abundance of the galaxy derived 
from our spectroscopic data, see Paper II) for all models. 

We found a relatively good correspondence (see Figs. 37, 
38 and 39, top, in Paper I) between the optical/A^ZR data 
and the models, especially for compact and dwarf objects 
such as HCG 31 Fl or SBS 0948-^532, the ages being con- 



sistent with a recent star-formation event (<100 Myr). We 
remark here 

1. the quality of the observational data and the data re- 
duction process, which was performed in detail and in 
an homogeneous way for all galaxies, 

2. the we corrected the data for extinction and reddening, 
considering the c(II/3) value derived from the spectro- 
scopic data obtained for each region (Paper II). As we 
have seen, this correction is important and very often it 
is not performed in the analysis of the colours of extra- 
galactic objects, which only consider the extinction of 
the Milky Way in the direction to the analysed galaxy, 

3. and the correction of the colours for the gas emission 
using our spectroscopic data. This effect is not impor- 
tant in some galaxies, but it seems fundamental when 
analysing compact objects with strong nebular emis- 
sion, such as BCDGs or regions within a galaxy pos- 
sessing an strong star burst. 

Some discordances between the colours and the predictions 
of the theoretical models ('^0.2 mag or even higher) are al- 
ways found in galaxies hosting a considerable population of 
old stars (Mkn 1199, Mkn 5, Tol 1457-262 #15 and #16, 
ESO 566-7), because their luminosities barely contribute to 
the U magnitude. Hence, the young stellar population usu- 
ally dominates the U — B colour, but the rest of the colours 
{B — V, V — R, V — J, H — Kg) posses an important contri- 
bution of the old stellar population. If the bursts and the 
underlying stellar population are analysed independently, 
the agreement between colours and the predictions given 
by the models is closer that when considering the galaxies 
as a whole. The last columns in Table [1] compile the ages 
estimated for the most recent star-forming event and the 
underlying population component (if possible) for all indi- 
vidual galaxies derived from optical/A7i? colours. Because 
the theoretical models are optimized to study the youngest 
stellar populations within the galaxies, in some cases we 
considered a lower limit of 500 Myr for the age of the un- 
derlying component (UC). 

Figure [5] shows several colour-colour diagrams compar- 
ing the predictions given by evolutionary synthesis models 
with the colours (corrected for both reddening and contri- 
bution of the emission lines) of our galaxy sample when 
the burst (blue stars) and underlying component (UC, red 
squares) of each system are considered independently. The 
correspondence has improved now. Indeed, all inferred ages 
of the most recent star-formation burst are lower than 
25 Myr, while the data corresponding to the underlying 
component suggest ages higher than 100 Myr. Therefore, a 
proper estimate of the stellar population age for this type 
of galaxy using broad-band filters is only obtained when 
bursts and underlying components are independently con- 
sidered. 

There are still some discrepancies that can be explained 
by a lack of good separation between regions with and with- 
out star-formation activity. The comparison of the V — R 
colour vs. the U — B colour in data of the UC also suggest 
that in some galaxies the old stellar population colours are 
not explained by just one single-age population, but at least 
two of them are needed (i.e., for SBS 1415+437, the UC 
colour may be explained by a mix of two stellar populations 
with ages of ~150 Myr and >500 Myr, see Sect. 3.15.1 in 
Paper I). However, the best method to analyse the colours 
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BC03 models. 



and luminosities of the host component in starburst sys- 
tems (specially, in BCGs) is performing a careful 2D anal- 
ysis of their structural parameters (i.e. Amorm et al. 2007; 
2009). Some of the galaxies analysed by these authors were 
also studied here. Their results of the colours of the UC (the 
host) agree well within the errors with those estimated here, 
for example, for Mkn 5 they compute (B-T^)t/c=0.40±0.28 
and {V — i?)t/c=0.28±0.16, while in this work we derived 
{B - T/)c/c=0.45±0.08 and {V - i?)c/c=0.30±0.08 for the 
same object. 

Figure [3] plots the Ha equivalent width (obtained from 
our narrow-band images) as a function of the U — B colour 
(obtained from our broad-band images) for the bursts 
within the analysed galaxies. We remark that the M^(Ha) 
derived from the Ha images agree quite well with those 
obtained from the optical spectroscopy (see Paper H). 
This figure compares t he observational data with some 
Starburst 99 models (jLeitherer et al.|[l999l ) at different 
metaUicities. As we see, the agreement is quite good for al- 
most all objects. This also indicates both the quality of our 
data and the success of the theoretical models to reproduce 
the young star-forming populations. 

However, if we compare the VF(Ha) and the U—B colour 
considering the total extension of each galaxy (and not only 
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the burst component), this agreement is less good. In this 
case, data with a fixed T/F(Ha) have a redder U — B colour 
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Table 2. Physical properties of the ionized gas for the galaxies analysed in this work. 



Galaxy 


Tc 


To High 


Tq Low 


Tie 


c(H/3) 


' * aos 


-W(K0) 




(n) 






[ciii J 






fAl 




HCG 31 AC 


D 


9400+600 


10800+300 


210+70 


0.09+0, 


.03 


2.0 


91.1+2.1 


HCG 31 B 


D 


11500+700 


12000+400 


<100 


0.28+0, 


.08 


2.0 


12.9+0.5 


HCG 31 E 


D 


11100+1000 


11800+600 


<100 


0.11+0, 


.05 


2.0 


21.1+1.1 


HCG 31 Fl 


D 


12600+1100 


12600+700 


<100 


0.32+0, 


.06 


2.0 


218+13 


HCG 31 F2 


D 


12300+1300 


12400+800 


<100 


0.14+0, 


,05 


2.0 


256+30 


HCG 31 G 


D 


11600+700 


12000+400 


<100 


0.09+0, 


.05 


2.0 


37.0+1.6 


Mkn 1087 


EC 


7100+1000 


8000+1000 


220+50 


0.17+0, 


.02 


1.7+0.2 


22.3+0.9 


Mkn 1087 N 


EC 


10900+1000 


10600+1000 


115+50 


0.17+0, 


,02 


0.2+0.1 


25.0+1.7 


Haro 15 C 


EC 


9500+800 


9600+600 


<100 


0.11+0, 


,03 


2.4+0.4 


16.4+1.1 


Haro 15 A 


D 


12850+700 


12000+500 


<100 


0.33+0, 


.03 


1.3+0.3 


75.7+4.2 


Mkn 1199 


D 


5400+700 


6800+600 


300+100 


0.30+0, 


,03 


1.8+0.4 


21.4+1.3 


Mkn 1199 NE 


EC 


8450+800 


8900+600 


<100 


0.16+0, 


,03 


0.6+0.3 


20.2+2.3 


Mkn 5 


D 


12450+650 


11700+450 


<100 


0.17+0, 


.02 


0.8+0.2 


75+5 


IRAS 08208+2816 


D 


10100+700 


10100+500 


<100 


0.11+0, 


,02 


3.2+0.1 


80+5 


IRAS 08339+6517 


EC 


8700+1000 


9100+1000 


100 


0.22+0, 


,02 


1.8+0.2 


19.0+0.8 


IRAS 08339+6517C 


EC 


9050+1000 


9300+1000 


100 


0.18+0, 


.03 


1.5+0.2 


7.5+0.2 


POX 4 


D 


14000+500 


12800+400 


250+80 


0.08+0, 


,01 


2.0+0.1 


200+9 


POX 4 Comp 


EC 


12400+1000 


11700+700 


<100 


0.06+0, 


,03 


1.4+0.3 


12+4 


UM 420 


D 


13200+600 


12200+500 


140+80 


0.09+0, 


.01 


2.0+0.1 


169+10 


SBS 0926+606 A 


D 


13600+700 


12500+500 


<100 


0.12+0, 


,03 


0.7+0.1 


125+6 


SBS 0926+606 B 


EC 


11500+1000 


11000+800 


<100 


0.18+0, 


,04 


1.0+0.3 


18+3 


SBS 0948+532 


D 


13100+600 


12200+400 


250+80 


0.35+0, 


.03 


0.3+0.1 


213+11 


SBS 1054+365 


D 


13700+900 


12600+700 


<100 


0.02+0, 


.02 


0.8+0.1 


89+7 


SBS 1054+365 b 


EC 


11800+1100 


11300+900 


300+200 


0.03+0, 


,03 


0.3+0.1 


8+3 


SBS 1211+540 


D 


17100+600 


15000+400 


320+50 


0.12+0, 


.01 


1.3+0.1 


135+10 


SBS 1319+579 A 


D 


13400+500 


12400+400 


<100 


0.03+0, 


.01 


0.0+0.1 


285+14 


SBS 1319+579 B 


D 


11900+800 


11300+600 


<100 


0.11+0, 


,03 


0.4+0.1 


42+4 


SBS 1319+579 C 


D 


11500+600 


11050+400 


<100 


0.02+0, 


.02 


0.2+0.1 


94+6 


SBS 1415+579 C 


D 


16400+600 


14500+400 


<100 


0.01+0, 


.02 


0.8+0.1 


222+11 


SBS 1415+579 A 


D 


15500+700 


13850+500 


<100 


0.15+0, 


,03 


1.0+0.2 


130+8 


III Zw 107 A 


D 


10900+900 


10500+800 


200+60 


0.68+0, 


.04 


2.0+0.3 


44+2 


Tol 9 


D 


7600+1000 


8300+700 


180+60 


0.50+0, 


.05 


7.5+0.8 


33+2 


Tol 1457-262 A 


D 


14000+700 


12500+600 


200+80 


0.57+0, 


,03 


1.4+0.2 


101+6 


Tol 1457-262 B 


D 


15200+900 


14200+700 


<100 


0.00+0, 


,05 


0.0+0.1 


82+7 


Tol 1457-262 C 


D 


13400+1100 


12400+1000 


200+100 


0.15+0, 


.02 


0.7+0.1 


92+9 


ESQ 566-8 


D 


8700+900 


9100+800 


300+100 


0.49+0, 


,03 


1.3+0.1 


95+7 


ESO 566-7 


EC 


7900+1000 


8500+900 


100+50 


0.23+0, 


,05 


2.7+0.2 


13+2 


NGC 5253 A 


D 


12100+260 


11170+520 


580+110 


0.23+0, 


.02 


1.3+0.1 


234+5 


NGC 5253 B 


D 


12030+260 


11250+490 


610+100 


0.38+0, 


,03 


1.7+0.1 


254+5 


NGC 5253 C 


D 


10810+230 


10530+470 


370+80 


0.25+0, 


,03 


0.8+0.1 


94+3 


NGC 5253 D 


D 


11160+510 


10350+650 


230+70 


0.10+0, 


.02 


0.6+0.1 


39+2 



° In this column we indicate if To was computed using the direct method (D) or via empirical calibrations (EC). 



than that predicted by the models. This is explained be- 
cause the U — B colour is slightly contaminated with the 
light of older stellar populations or regions with no nebular 
emission. Hence, as we emphasized before, it is important to 
distinguish between the pure starburst regions and the un- 
derlying component to get a good estimate of the properties 
of these galaxies and, in particular, the strong star-forming 
regions. 

The age of the last starburst event experienced by each 
galaxy and the minimum age of its old stellar popula- 
tions are compiled in the last two columns of Table [T] 
Except for a few objects (HCG 31 members E, Fl and 
F2 and Mkn 1087 members N and #1) for which it was 
not possible to estimate the colours of the UC, all anal- 
ysed galaxies show an older stellar population underly- 
ing the bursts. Indeed, in many cases the colours of the 
UC suggest ages older than 500 Myr. This clearly indi- 
cates that all galaxies have experienced a previous star- 
formation events long time before those they are now host- 
ing, as concluded in many other previous results (i.e. Cairos 
et al. 2001a,b; BergvaU & Osthn, 2002; Papaderos et al. 



2006; Amorm et al. 2009). However, as we pr eviously said 
(jLopez-Sanchez. Esteban k, Rodn'gue9 l2004al) . this seems 
not to be true in the particular case of members F of 
HCG 31, which clearly show no evidences of underlying old 
stellar populations. This was rece ntly confirmed by dee p 
Hubble Space Telescope imaging (iGallag her et all boiOh . 
and hence these two objects are very likely experiencing 
their very first star-formation event. 

3. Physical properties of the ionized gas 

Table [2] compiles all the high- and low-ionization electron 
temperatures of the ionized gas, T^, electron density rig, red- 
denning coefficient c(H/3) , equivalent width of the underly- 
ing stellar absorption in the Balmer H i lines Wabs, and the 
H/3 equivalent width, for the galaxies analysed in this work 
(see Paper II). Thirty-one up to 41 of the objects listed in 
Table [5] have a direct estimate of the electron temperature 
of the ionized gas. For most, this was computed using the 
[Oiii] ratio involving the nebular [Oiii] AA4959,5007 and 
the auroral [O ill] A4363 emission lines. In most objects, 
the low-ionization electron temperature was not computed 
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Fig. 4. Comparison of some observational flux ratios obtained 
for all available regions analysed in this work with the diagnostic 
diagrams proposed by Dopita et al. (2000), blue continuous line 
(DOO), and Kewley et al. (2001), red discontinuous line (KOI). 
The left panel also shows the empirical relation provided by 
Kauffmann et al. (2003) with a dotted-dashed dark yellow line 
(Ka03). 



directly but assumi ng; the re l ation between re(Oiii) and 
Te{0 II) provided bv lGarnettI (|1992[) . Half of the objects of 
our galaxy sample (22) have electron densities lower than 
100 cm-3. 

We explored possible correlations among some of the 
different quantities compiled in Table [5] as well as the oxy- 
gen abundance (see Table [3]) computed for the objects. 
First, we checked the nature of the ionized gas of the 
sample galaxies. Figure H] plots the typical diagnostic di- 
agrams between bright emission lines and t he predictions 
given by the photoionized models provided bv lDopita et al.l 
(j2000t) for extragalactic H ii regions (that assume instanta- 
neous st ar-formation within star-forming regions) and the 
iKewlev et al. (2001.) models for starburst galaxies (which 
consider continuous star formation and more realistic as- 
sumptions about the ph ysics of starburst ga laxies) . The di- 
viding line given by the lKewlev et al.l ()2001[ ) models repre- 
sents an upper envelope of positions of star-forming galax- 
ies. As we see, in all cases the data are found below the 
theoretical prediction given by this line. This indicates that 
photoionization is the main excitation mechanism of the 
gas. We will get the same result when we compare the re- 
lation between the FIR and the radio-continuum luminosi- 
ties (Paper V). It is interesting to notice that the obser- 
vational points included in the diagnostic diagram involv- 
ing the [Oiii]/H/? and [Nil] /H a ratios are l ocated close 
to the prediction given by the iDooita et all (|2000t ) mod- 
els, while points included in the diagnostic diagram that 
considers the [Oiii]/H/3 and [S ii] /Ha ratios are foti r id ver y 
close to the upper envelope given bv IKewlev etaU (|200ll) . 
The left panel in Figure |4] includes the empirical relation 
be tween the [Olll]/H/3 and the [NiiJ/Ha ratios provided 
bv iKauffmann et al.l (|2003l ) analysing a large data sample 
of star-forming galaxies from the Sloan Digital Sky Survey 
(SDSS; Y ork et al. 2000). As can be seen, the compari- 
son of the IKauffmann et al] (|2003f ) relation with our data 
points also indicates that our objects are experiencing a 
pure star- formation event, despite a clear offset between 
both datasets. 
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Fig. 5. High ionization electron temperature, re(Olll), vs. the 
H/3 equivalent width. The colour-dots indicate the oxygen abun- 
dance range of each object. Some objects have been labeled. 



W{H/3) is a good indicator for the age of the most re- 
cent star-formation event. The hydrogen ionizing flux of a 
star cluster gradually decreases as the most massive stars 
disappear with time, and hence the width of H/3 decreases 
with time (see Papers I, H and HI). Figure [5] plots the rela- 
tion between the high ionization electron temperature and 
the H/3 equivalent width. As we see, |W^(H/3)| increases with 
increasing Tq, but we must remember that there is a strong 
correlation between the electron temperature and the oxy- 
gen abundance, as high-metallicity H ii regions cool more 
efficiently than low-metallicity H ii regions. To study this 
effect, we used colours to plot four metallicity ranges in 
Fig. 121 These colours indicate the oxygen abundance range 
of each object: <7.90, 7.90-8.15, 8.15-8.40 and >8.40, in 
units of 12-|-log(0/H). Although now it is not so evident, 
it still seems that regions with larger |M^(H/3)| tend to have 
higher Tc- This indicates that younger bursts have a larger 
ionization budgets and arc therefore capable to heat the 
ionized gas to higher electron temperatures. Another ef- 
fect that we should considered here is that galaxies with 
higher metallicity (and hence with lower electron tempera- 
ture) usually have a higher absorption in the H/3 line than 
low-metallicity objects because of older stellar populations. 
Indeed, it is interesting to note that objects in the metallic- 
ity range 7.90<12-|-log(O/H)<8.15 may have any value of 
W{Il/3). This is very probably because within this metal- 
licity range lie both dwarf objects with no underlying old 
stellar population (i.e., HCG 31 F, SBS 0948-^532) and 
galaxies which possess a considerable amount of old stars 
(i.e, HCG 31 B, Mkn 5). 

Figure ini shows that the ionization degree of the ionized 
gas of the starburts seems to increase with the increas- 
ing of |iy(H/3)|. This is a relation similar to that found 
in Fig. IHI and also indicates that younger bursts harbour a 
higher proportion of massive stars and therefore their asso- 
ciated H II regions have larger ionization parameters. That 
is evident in NGC 5253 A, B and HCG 31 F, which show 
both the highest values of the |Ty(H/3)| and the 0++/0+ 
ratio and possess the youngest star-formation bursts (see 
Table [T|). The increasing of the 0"'""'"/0"'" ratio as increas- 
ing I (H/3) I seems to be independent of the metallicity. 
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Fig. 6. 0++/0+ ratio vs. W(H/3). The colour-dots indicate the 
oxygen abundance range of each object. Some objects have been 
labeled. 
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Fig. 7. Reddening coefficient, c(H/3), vs. the oxygen abundance 
for the regions analysed in this work. Some objects have been 
labeled. 



although galaxies with higher metallicity tend to show the 
lowest ionization degrees. 

The dependence of the reddening coefficient as a func- 
tion of other parameters is also interesting. Figure [7] plots 
c(H/3) vs. the oxygen abundance. Although the dispersion 
of the data is rather scattered, we see a clear dependence: 
the reddening coefficient is higher at higher metallicities. 
We should expect this result, because galaxies with higher 
oxygen abundance are chemically more evolved and should 
contain a larger proportion of dust particles that absorbs 
the nebular emission. 

Another interesting relation is shown in Fig. [SJ which 
plots the equivalent width of the stellar absorption un- 
derlying the Hi Balmer lines (Wabs) as a function of the 
oxygen abundance. We can see that objects with higher 
metallicities show larger Wabs- More metallic galaxies cor- 
respond to more massive and chemically evolved systems, 
which means that they have consumed a larger fraction of 
their gas and the stellar component should be compara- 
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Fig. 8. Equivalent width of the underlying stellar absorption 
in the Balmer Hi lines, Wats, vs. the oxygen abundance for 
the regions analysed in this work. Notice that the y-axis has 
been broken for clarity. The value of Wats for UM 420 is very 
probably overestimated, because this galaxy is observed through 
the external areas of the spiral disc of the foreground galaxy 
UCG 01809 (see Fig. 13 in Paper I). 



tively more important. The data corresponding to the low- 
est metallicity objects analysed in this work (SBS 1415-1-579 
and SBS 1211-1-540) show a value of Wabs relatively high 
what is to be for them. This suggests a considerable under- 
lying stellar population in these very low-metallicity galax- 
ies, as we already discussed (see Sect. 3.15 and Sect 3.13 of 
Paper I). 

Figure [9] plots VF(H/3) vs. the oxygen abundance. The 
very large dispersion of W{Hf5) for 12-t-log(0/H) of about 
8.0 is remarkable, but it also seems clear that galaxies with 
0/H ratios higher than that value tend to have |VF(H/3)| < 
100 A and, conversely, galaxies with lower oxygen abun- 
dances show |Ty(H/3)| > 100 A. This behaviour may be 
be related to the results in Fig. |8l in the sense that more 
metallic objects tend to have a higher underlying stellar ab- 
sorption that can lead to an underestimation of |iy(H/3)|. 
Indeed, Fig. |9] also compares our observational data with 
the predictions given by th e chemical evolution model s 
of Hii galaxies provided bv iMartm-Manion et al.l (|2008[ ). 
They assumed the star formation as a set of successive 
bursts, each galaxy experiencing 11 star-formation bursts 
along its evolution of 13.2 Gyr. Figure M includes the re- 
sults for the first {t—0 Gyr), second (t=1.3 Gyr) and last 
(i=13.2 Gyr) bursts for a model that considers an atten- 
uated bursting star- formation mode and that 1/3 of the 
gas is always used to form stars in each time-step. As we 
see, all the strong starbursting systems are located between 
the positions of the first and second burst models, con- 
firming that although the dominant stellar population is 
certainly very young, previous star-formation events in the 
last 500-1000 Myr are needed to explain our observational 
data points. This agrees well with the minimum ages of the 
underlying stellar component we derived using our photo- 
met ric_data^_(see_Sect. 2 and last column of Table [1]). The 
M artm-Manion et al.l (j2008f) models also explain the large 
dispersion of W(Rl3) for 12-hlog(0/H) of about 8.0, as well 
as the trend that more metal rich galaxies have lower values 
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Fig. 9. H/3 equivalent width vs. oxygen abundance. Some ob- 
jects have been label ed. Some chemical evo l ution models of H ll 
galaxies provided bv lMartm-Manion et alJ l|2008f ) are also plot- 
ted. See text for details. 
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Fig. 10. F([Oni] A5007)/F(H^) vs. W(H^) and comparison 
with the predictions given by four different ionization models 
bv fStasihska et al.. (,2001; '). See text for details. 



of |T^(H/3)| because of the effects of the underlying stellar 
populations. 

Another indication of the effect of the underlying 
evolved stellar population is found in Fig. [TUl which 
compares the [Olll] A5007 line flux with the W(Iip) 
of the sample galaxies with the rn odel predictions by 
IStasihska. Schaerer fc Leithereil (|200lh . Although a general 
good correspondence is found, some of the objects are 
slightly displaced to the left -lower W{B.p)- of the models 
predictions, suggesting that perhaps the measured values 
of |M^(H/3)| are underestimated for some of them that pre- 
cisely coincide with those with a larger oxygen abundance, 
as we also concluded before. 



4. Chemical abundances of the ionized gas 

Table [3] compiles the oxygen abundance, the 0++/0+ ratio 
and the N/0, S/0, Ne/0, Ar/0, and Fe/0 ratios (ah in 
logarithmic units) for all galaxies analysed in this work. As 
we already said, 31 up to the 41 independent regions within 
the WR galaxy sample analysed here have a direct esti- 
mate of the electron temperature (as indicated in Table [31). 
Figure [TT] (left) plots the N/0 ratio vs. the oxygen abun- 
dance for all our data with a direct estimate of the electron 
temperature and its comparison wit h previous samples in- 
volving similar obj ects and To-based (jlzotov fc Thuaniri999l : 
llzotov et ani2004li . This figure shows that the position of 
our data agrees with that obtained using other observa- 
tions. The errors we estimated in our objects are in gen- 
eral higher than t hose r eported by llzotov fc ThuanI (jl999t) 
and llzotov et al.l ()2004D basically because we used differ- 
ent criteria for estimating observational errors, which are 
more conservative as well as more realistic in our opin- 
ion. We also point out that our data are always of higher 
spectral and spatial resolution than those obtained by the 
aforementioned authors, and have a similar or even higher 
signal-to-noise ratio in many cases. This is an important 
point to be clarified because non-specialist in the spectra 
of ionized nebulae may interpret that lower quoted uncer- 
tainties are synonymous of better observational data, and 
this may not always be the case. Although the errors in the 
electron temperatures derived using the empirical methods 
are large, relative atomic abundances (such as the N/0 ra- 
tio) are less sensitive to the choice of To- Therefore they 
are used in many occasions to compare with the results 
provided by previous observations or with the predictions 
given by theoretical models. 

4.1. The N/0 ratio 

Figure |TT] (right) plots the N/0 ratio vs. 12-Klog(0/H) 
for all objects analysed in this work; the chemical abun- 
dances were derived either from the direct method or 
via empirical calibrations. We compare o ur data with the 
two galaxy samples previously indicated (jlzotov fc Thuai] 
1 19991 : llzotov et all |2004[ ) and with other galaxy samples 
whos e data have been o btained using empirical calibra- 
tions: llzotov et al.l (l2006^■ whose data were ext racted from 
the SPSS, and lvan Zee. Salzer fc Havnei ()1998[ ). who study 
data from H ii regions within spiral galaxies with chem- 
ical abu ndances co i npute d via the direct method or us- 
ing the iMcGaughl ()1994[ ) empirical calibration. In some 
sense, the N/0 ratio of a galaxy is an indicator of the 
time that has elapsed since the bulk of star formation oc- 
cu rred, or of the nominal a ge of the galaxy as suggested 
bv lEdmunds fc Paeell (|1978[ ). Following the position of our 
data points in Fig. [TT]we see that they follow the expected 
trend: 

1. The N/0 ratio is rather constant for 12+log(0/H)<7.6. 
In our case, for the galaxies SBS 1211-1-540 and 
SBS 1415-1-437 we deri ve log N/0~ -1 . 6, sim ilar val- 
ues as those found by llzotov fc ThuanI (|1999^ . These 
authors explained the constant N/0 ratio in very low- 
metallicity objects assuming that the nitrogen is pro- 
duced only as a primary element in massive, short- 
life stars. However, other authors have claimed that 
this may be not completely true (i.e., Henry et al. 
2000; Pilyugin et al. 2003; Molla et al. 2006) because 
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Table 3. Chemical properties of the ionized gas for the galaxies analysed in this work. 



Galaxy Tc" 12+log(0/H) log i^-q- log(N/0) log(S/0) log(Nc/0) log{Ar/0) log{Fe/0) 



HCG 31 AC 


D 


8 


22+0.05 


1.51+0.12 


-1 


12+0 


08 








-0 


93+0.12 




-2 


12+0 


21 


HCG 31 B 


D 


8 


14+0.08 


0.63+0.09 


-1 


39+0 


10 


-1 


67+0 


14 


-0 


42+0.13 




-1 


87+0 


32 


HCG 31 E 


D 


8 


13+0.09 


1.00+0.11 


-1 


26+0 


12 


-1 


58+0 


15 


-0 


42+0.14 




-1 


77+0 


32 


HCG 31 Fl 


D 


8 


07+0.06 


3.72+0.32 


-1 


27+0 


11 


-1 


69+0 


15 


-0 


80+0.17 






-1.9: 




HCG 31 F2 


D 


8 


03+0.10 


2.19+0.21 


-1 


43+0 


16 


-1 


67+0 


18 


-0 


76+0.20 






... 




HCG 31 G 


D 


8 


15+0.07 


1.15+0.11 


-1 


31+0 


10 


-1 


67+0 


22 


-0 


56+0.14 






-2.0: 




Mkn 1087 


EC 


8 


57+0.10 


0.55+0.18 


-0 


81+0 


12 


-1 


78+0 


16 


-0 


45+0.17 










Mkn 1087 N 


EC 


8 


23+0.10 


0.99+0.25 


-1 


46+0 


15 








-0 


52+0.19 










Haro 15 C 


EC 


8 


37+0.10 


-0.23+0.16 


-1 


03+0 


15 


-1 


71+0 


18 


-0 


65+0.18 






-2.2: 




Haro 15 A 


D 


8 


10+0.06 


0.66+0.10 


-1 


35+0 


11 


-1 


89+0 


15 


-0 


68+0.12 






-1.6: 




Mkn 1199 


D 


8 


75+0.12 


-0.36+0.16 


-0 


62+0 


10 


-1 


54+0 


14 


-0 


58+0.17 




-1 


86+0 


26 


Mkn 1199 NE 


EC 


8 


46+0.13 


-0.19+0.09 


-1 


20+0 


11 


-1 


54+0 


17 


-0 


65+0.18 










Mkn 5 


D 


8 


07+0.04 


0.25+0.08 


-1 


38+0 


07 


-1 


62+0 


11 


-0 


80+0.13 


-2.31+0.12 


-1 


96+0 


18 


IRAS 08208+2816 


D 


8 


33+0.08 


0.43+0.12 


-0 


89+0 


11 


-1 


64+0 


16 


-0 


67+0.13 


-2.51+0.15 


-1 


95+0 


17 


IRAS 08339+6517 


EC 


8 


45+0.10 


0.53+0.16 


-0 


94+0 


14 








-0 


45+0.18 










IRAS 08339+6517C 


EC 


8 


38+0.10 


0.81+0.21 


-1 


13+0 


17 










-0.55: 










POX 4 


D 


8 


03+0.04 


0.74+0.06 


-1 


54+0 


06 


-1 


80+0 


10 


-0 


78+0.10 




-2 


17+0 


11 


POX 4c 


EC 


8 


03+0.14 


-0.30+0.22 


-1 


60+0 


20 










-0.60: 










UM 420 


D 


7 


95+0.05 


0.00+0.08 


-1 


11+0 


07 


-1 


66+0 


13 


-0 


71+0.13 




-2 


16+0 


13 


SBS 0926+606 A 


D 


7 


94+0.08 


0.42+0.12 


-1 


45+0 


09 


-1 


60+0 


13 






-2.34+0.13 


-1 


99+0 


16 


SBS 0926+606 B 


EC 


8 


15+0.16 


0.21+0.14 


-1 


35+0 


12 




















SBS 0948+532 


D 


8 


03+0.05 


0.61+0.08 


-1 


42+0 


08 


-1 


69+0 


14 


-0 


73+0.12 




-1 


78+0 


10 


SBS 1054+365 


D 


8 


00+0.07 


0.70+0.11 


-1 


41+0 


08 


-1 


79+0 


15 


-0 


67+0.11 


-2.29+0.14 








SBS 1054+365 b 


EC 


8 


13+0.16 


-0.35+0.20 


-1 


47+0 


20 




















SBS 1211+540 


D 


7 


65+0.04 


0.69+0.07 


-1 


62+0 


10 


-1 


47+0 


14 


-0 


75+0.10 










SBS 1319+579 A 


D 


8 


05+0.06 


0.77+0.12 


-1 


53+0 


10 


-1 


76+0 


10 






-2.41+0.11 








SBS 1319+579 B 


D 


8 


12+0.10 


0.16+0.19 


-1 


49+0 


12 


-1 


76+0 


14 














SBS 1319+579 C 


D 


8 


15+0.07 


0.18+0.13 


-1 


38+0 


10 


-1 


60+0 


11 










-2.3: 




SBS 1415+437 C 


D 


7 


58+0.05 


0.35+0.08 


-1 


57+0 


08 


-1 


62+0 


12 






-2.31+0.13 


-1 


91+0 


13 


SBS 1415+437 A 


D 


7 


61+0.06 


0.42+0.14 


-1 


57+0 


09 


-1 


72+0 


14 










-1.9: 




ITT Zw 107 A 




g 




n 1 2+0 1 4 


_x 


16+0 


10 


_l 


82+0 


15 


-0 


73+0.15 


_o 4f;+n 1 3 




-2.3: 




Tol 9 


D 


8 


57+0.10 


0.16+0.17 


-0 


81+0 


11 


-1 


62+0 


12 


-0 


72+0.14 


-2.55+0.15 




-2.1: 




Tol 1457-262 A 


D 


8 


05+0.07 


0.27+0.11 


-1 


57+0 


11 


-1 


88+0 


13 


-0 


88+0.18 


-2.50+0.13 




-2.2: 




Tol 1457-262 B 


D 


7 


88+0.07 


0.43+0.11 


-1 


61+0 


12 


-1 


72+0 


18 


-0 


88+0.20 


-2.44+0.18 


-1 


90+0 


22 


Tol 1457-262 C 


D 


8 


06+0.11 


0.14+0.16 


-1 


59+0 


16 








-0 


84+0.22 


-2.45+0.20 








ESQ 566-8 


D 


8 


46+0.11 


-0.19+0.17 


-0 


76+0 


12 








-0 


56+0.19 


-2.17+0.19 




-2.5: 




ESQ 566-7 


EC 


8 


50+0.16 


-0.57+0.22 


-0 


82+0 


16 












-2.49+0.25 








NGC 5253 A 


D 


8 


18+0.04 


2.88+0.18 


-0 


91+0 


07 


-1 


58+0 


08 


-0 


71+0.08 


-2.19+0.07 


-2 


10+0 


12 


NGC 5253 B 


D 


8 


19+0.04 


3.09+0.14 


-1 


02+0 


07 


-1 


60+0 


08 


-0 


70+0.08 


-2.21+0.07 


-2 


18+0 


11 


NGC 5253 C 


D 


8 


28+0.04 


1.95+0.13 


-1 


50+0 


08 


-1 


69+0 


09 


-0 


74+0.08 


-2.30+0.08 


-2 


46+0 


14 


NGC 5253 D 


D 


8 


31+0.07 


0.56+0.14 


-1 


49+0 


10 


-1 


74+0 


13 


-0 


70+0.15 


-2.30+0.13 


-2 


25+0 


16 



" In this column we indicate if To was computed using the direct method (D) or via empirical calibrations (EC). 



of the lack of a clear mechanism that produces N in 
massive stars besides the effect or the stellar rotation 
([Mevnet fc Macdcr 2005V Furthermore, these galaxies 
already host old stellar populations, and hence low- and 
intermed iate-mass stars sho uld be also releasing N to 
the ISM. lHenrv et all (12000( 1 explained the constancy of 
the N/0 ratio in metal-poor galaxies by a historically 
low star-formation rate, where almost all the nitrogen 
is pro duced by 4-8 Mq stars. Additionally. Ilzotov et al.l 
(|200(3[ ) suggested that the low dispersion of the data in 
this part of the diagram is probably explained by the 
low number of WR stars that are expected at very low- 
metallicity regimes. 
2. However, there is an important dispersion of the data 
in the interval 7.6<12+log(0/H)<8.3. This problem 
has been analysed by several authors in the past (i.e., 
Kobulnicky & Skillman, 1998; Izotov & Thuan, 1999; 
Pilyugin et al. 2003; Molla et al. 2006). Two main sce- 
narios have been proposed for explaining this dispersion: 
(a) A loss of heavy element via galactic winds. In partic- 
ular, it should be a loss of a-elements via supernova 
explosions, a-elements, such as oxygen, are produced 



in massive short-lived stars ()Edmunds fc Pagelll978l : 
IClavton fc Pantelakilll993l ). Hence, the effect of su- 
pernova explosion s would produce a underabun- 
dance of oxygen (|Esteban fc PeimbertI ll995), in- 
creasing the N/0 ratio. However, the observational 
evidence for low-mass galaxies with galactic winds 
is still unclear (i.e. Marlowe et al. 1995; Bomans et 
al. 2007; Dubois fc Teyssier 2008; van Eymeren et 
al. 2008, 2009, 2010) and even the numerical simu- 
lations give very discrepant results (i.e., Mac Low & 
Ferrara 1999; Springel & Hernquist 2003; Tenorio- 
Tagle et al. 2006; Dubois & R. Teyssier 2008). 
(b) A delayed release of nitrogen and elements pro- 
duced in low-mass long-lived stars compared to the 
a-elements. The N/0 ratio drops and the 0/H ratio 
increases as supernovae release the a-elements into 
the ISM. Consequently, the chemical properties of 
these galaxies would vary ve ry fast (few tens of Myr) 
durin g the starburst phase (jKobulnickv fc SkillmanI 
I1998D . The delayed- release hypothesis also predicts 
that BCDGs with high N/0 ratios are experienc- 
ing their first burst of massive star formation after 
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Fig. 11. [Left) N/0 ratio vs. oxygen abundance for all objects with a direct estimate of Tc. Our data are compared with those 
compiled by IT99 (Izotov & Thuan 1999) and Izotov et al. (2004). [Right) N/0 ratio vs. oxygen abundance for all galaxies analysed 
in this work (also included objects which chemical abundances are derived using empirical calibrations). For comparison, we have 
included other galaxy samples from the literature whose O/H ratio has been determined using empirical calibrations: Izotov et al. 
(2006), van Zee et al. (1998), IT99 (Izotov & Thuan 1999), Izotov et al. (2004) and PTV03 (Pilygin et al. 2003). 



a relatively long quiescent interval (oxygen has still 
not been completely delayed by massive stars and 
mixed with the surrounding ISM), while those ob- 
jects with low N/0 ratios have had little or no qui- 
escent interval. However, recent chemical evolution 
models suggest that the releasing and mixing of the 
oxygen occurs almost instantaneously, and hence the 
delayed-release scenario cannot explain BCDGs with 
high N/0 ratio. If this is the case, the most plausible 
explanation of the high N/0 ratio observed in these 
objects is the chemical pollution due to the winds of 
WR stars, which are indeed ejecting N to the ISM, 
as we will discuss below. 
3. For moderate high-metallicities objects, 
12-|-log(0/H)>8.3, the N/0 ratio clearly increases 
with increasing oxygen abundance. This trend seems 
to be a consequence of the metallicity-dependence of 
nitrogen production in both massive and intermediate- 
mass stars (e.g., Pilyugin et al. 2003), the N/0 ratio 
increases at higher metallicities. Hence, nitrogen 
is essentially a secondary element in this metal- 
licity regime (Torres-Pcimbert. Pei mbert fc Fierra 
1981 IVila-Costas fc Edmundsi 119931 



Henrv. Edmunds fc KoppenT 20001 Ivan Zee fc HavnesI 
2006t ). Besides the uncertainties (that are higher than 



those estimated in other objects because the Te error 
is higher at higher metallicities) our data agree with 
the tendency found in other galaxy samples, as that 
compiled by iPilvugin et al.l ( 20031) . H owever, notice 
that the galaxy sample compiled by Ivan Zee et al.l 
(|1998[ 1 does not agree with our data, as her data 
have a systematically lower N/0 ratio. In some cases, 
differences higher than 0.5 dex in the N/0 ratio 
are found for a particular oxygen abundance. This 
discrepancy may be partially explained by the fact 
t hat empirical calibr ations from phot o ionize d models 
- Ivan Zee etaD llMl) used iMcGauehl (|l99ll) models- 
seem to overestimate the actual oxygen abundance by 
at least 0.2 dex (see below), and hence the derived N/0 
ratio is lower than the actual value. 



In summary. Fig. [TT] can be explained assuming the very dif- 
ferent star-format ion histories that each individual galaxy 
has e xperienced (iPilvugin et al.l l2003t Ivan Zee fc HavnesI 
l2006t ). A galaxy with a constant SFR will have a lower net 
N/0 yield than a galaxy with declining SFR, because more 
oxygen has been released into the ISM due to the ongo- 
ing star-formation activity. This observational result com- 
pletely agrees with the predictions given by chemical evo- 
lution models that consider the effect of the star-formation 
history in the N/ O-O/H diagram, as those presented by 
iMoUa et al.l (pOOQ) . 



4.2. Nitrogen enrichment in WR galaxies 

From Fig. |Tl] it is evident that there are some objects 
in the low-intermediate metallicity regime with a higher 
N/0 ratio than expected for their oxygen abundance. An 
excess of nitrogen abundance has been reported in a few 
cases (e.g. Kobulnicky et al. 1997, Pustilnik et al. 2004). 
Remarkably, the common factor observed in all galaxies 
with a high N/0 ratio is the detection of Wolf-Rayet 
features. In deed, as we demonstrated in our analysis of 
NGC 5253 (jLopez-Sanchez et al.l[2007l ). the ejecta of WR 
stars may be the origin of a localized N (and probably also 
He) enrichment of the ISM in these galaxies. 

The analysis of the WR populations within our sam- 
ple galaxy was performed in Paper HI. The numbers of 
WNL and WCE stars were coi nputed assuming metallicity - 
dependent WR luminosities (jCrowther fc Hadfieldl 120061 ) . 
We detected the blue WR bump (the broad Hell A4686 
emission line) in all objects with high N/0 ratio: UM 420, 
NGC 5253 A,B, HCG 31 AC, IRAS 08208-f 2816, HI Zw 107 
and ESO 566-8, indicating that these bursts host an impor- 
tant WNL star population. 

Figure [l2| plots the observed nitrogen overabundance, 

A[N/0) = l0g(7V/0obscrvcd) - log(iV/Oexpcctod ft'om O/h), 

as a function of the derived WNL/(WNL-|-0) ratio. We do 
not see any clear trend in this diagram, merely that the 
objects with a high N/0 ratio do not show a particularly 
high WNL/(WNL-|-0) ratio. 
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WNL/(WNL+0) ratio for our sample galaxies (derived in 
Paper III). Some important objects have been labeled. 



Galaxies HCG 31 AC and III Zw 107 seem to show a 
slight nitrogen excess [A{N/0) ^ 0.15 de xj. Three of the 
;alaxi es with high N/0 ratio compiled by iPustilnik et all 
2004 are plotted in Fig.Hl NGC 5253 (already discussed), 
Mkn 1089 (HCG 31 AC) and UM 420. Our data confirm the 
nitrogen overabundance in UM 420 [A{N/0) ^ 0.4 dex], 
but not a significant N/0 ratio in HCG 31 AC [A{N/0) - 
0.15 dex: IPustilnik et alj (|2Q04f ) quoted - 0.5 dex]. Wc also 
find a relatively high N/0 ratio in ESO 566-8 [A{N/0) - 
0.3 dex]. 

The very rare occurrence of objects with a large N 
overabundance suggests the general idea of the short- 
time sc ales for the localized N pollut io n and its fast dis- 
persal. iBrinchmann. Kunth fc Durretj (|2008[ ) used SDSS 
data to find that for |1F(H/3)| > 100 A, WR galax- 
ies show a high N/0 compared to non-WR galax- 
ies. Quantitatively these authors found that on average 
Alog(N/O)[wR-nonWR]=0-13±0.04. They interpreted this 
result as a rapid enrichment of the ISM from WR winds. 
iBrin chmann ct al. (2008) also found that WR galaxies are 
in general more metal-rich at a given |Vl^(H/3)| than similar 
galaxies without WR features, which likely is a reflection of 
WR stars being more abundant at higher metallicities (see 
Fig. 5 of Paper HI). 

Althou gh we do not dismis s the s tatistical analysis per- 
formed by Brinchmann et aP (|2008[ ) comparing WR and 
non-WR galaxy data, we would like to warn about the use 
of data with low spectral resolution in order to derive an 
accurate nitrogen abundance in individual objects. This is 
commonly done via the analysis of the [N ii] A5683 emis- 
sion line, very close to Ha. Lack of sufficient spectral res- 
olution will derive a blending of both lines and a proba- 
ble over-estimation of the [N ii] A5683 flux, that would be 
even higher if broad low-intensity wings in the Ha pro- 
file exist, which are actually rather common in WR galax- 
ies (e.g. Men dez fc Esteban 199 7). On the other hand, as 
explained bv llzotov et al.l (|2006[ ). the bright doublet [On] 
AA3726,3729 is not observed in the SDSS spectra for nearby 
galaxies, and hence the estimate of the total oxygen abun- 
dance has to be done via the [O ii] AA7319,7330, which 
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Fig. 13. H/3 equivalent width vs. the N/O ratio. Blue circles 
indicate regions with no detected WR features. Red squares 
are galaxies in which WR features are observed. Yellow dia- 
monds correspond to objects with a high N/O ratio for their 
oxygen abundance, all showing WR features. The A{N/0) = 
log(Af/Oobservcd)-log(iV/Ocxpectod from o/h) difference is shown 
in every galaxy with dotted yellow lines. 



is much fainter, very dependent on the electron density 
and severely affected by sky emission features. Therefore 
their associated errors are usually larger than for the [O ii] 
AA3726,3729 lines. More data and a re-analysis of the chem- 
ical abundances in galaxies where WR features are de- 
tected, with a similar analysis of a sample of non-WR galax- 
ies, are needed to get any definitive results. 

Some authors (i.e. Izotov et al. 2006) suggest that 
there is a dependence between the N/O ratio and the 
H/3 equivalent width: the N/O ratio should increase 
wit h decreasing |iy(H / 3)|. T his trend was also observed 
by IBrinchmann et al] ()2008l ). but only for objects with 



|Vl^(H/3)| > 100 A. Figure[T3|plots the N/O ratio vs. iy(H/3) 
for the objects analysed here. In this figure we distinguish 
between objects with high N/O ratio for their oxygen abun- 
dance and WR features (yellow diamonds), and galaxies 
with a normal N/O ratio with (red squares) or without 
(blue circles) WR features. 

From Fig. \W\ it is evident that non-WR galaxies only 
show high N/O ratios when their |I^(H/3)| < 50 A. 
Their N /O ratio b ecomes low for larger equivalent widths. 
Brinch mann et al . (2008) suggested that the non-detection 
of high N/O ratios in objects with small equivalent widths 
is consistent with very young bursts where the WR stars 
have not yet had a change to enrich the surrounding ISM 
to a noticeable degree. However, this is probably a conse- 
quence of both the complex star-formation histories and 
the high relative importance of the old underlying stellar 
populations in these systems. 

On the other hand it is remarkable that the WR galax- 
ies with |W^(H;3)| <50 A show systematically high N/O 
ratios, but a large dispersion when |M^(H/3)| > 50 A. This 
dispersion becomes substantially smaller when we consider 
the N/O ratio expected for their O abundance for those 
objects with a high A(N/0). If the high N/O ratios in 
these galaxies are produced by the chemical pollution due 
to winds of WR stars, objects would move to the right in 
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Fig. 14. Ne/O, S/O, Ar/O ang Fe/O ratios vs the oxyg en abundance for the g alaxi es analysed in this work. We compare with 
the results provided by the galaxy samples considered bv llzotov fc ThuanI (|1999f ) and llzotov et al.l (|2006l ). The dashed dark- blue 
lines indicate a fit to our data. 



Fig. [T3] once the burst is finished because of the decreasing 
of |W(H/3)|. If the chemical poUution is very locahzed, we 
should also expect that objects with a N excess would move 
towards lower values of the N/0 ratio as the fresh released 
material is dispersed and mixed with the existing gas of 
the galaxy. Hence, detailed analyses of galaxies with a lo- 
calized high N/0 ratio, suc h as we performed in NGC 5253 
(jLopez-Sanchez et al.ll2007l ). are fundamental to solve these 
unsolved questions, which will definitely be key elements to 
the evolution of the galaxies. 

4.3. The a-elements to oxygen ratio 

Figure [H plots the Ne/O, S/O, Ar/O, and Fe/O ra- 
tios as a function of 12-|-log(0/H) for all objects anal- 
ysed in this work. We compare our observational results 
with those found in sarnples with similar characteris- 
tics (jlzotov fc Thuanlll999trizotov et al.ll2006D . We remark 
again that although the estimates of our errors are higher 
than those provided by other samples, they are not a con- 
sequence of the quality of our data but derive from the 
different formalism we used to estimate the uncertainties. 
In any case, our results completely agree with those pre- 
viously reported in the literature. The relative abundance 
ratios of the a-elements (neon, sulfur, and argon) to oxygen 
are approximately constant, as expected because all four 
elements (O, Ne, S and Ar) are mainly produced by mas- 



sive stars. For our sample, the mean log(Ne/0), log(S/0), 
and log(Ar/0) ratios are -0.70 ± 0.13, -1.68 ± 0.10 and 
—2.37 ± 0.12, respectively. These values are comparable to 
the values reported for starbursting dwarf galaxies (e.g., 
Izotov & Thuan, 1999; Izotov et al. 2006) and for other 
dwarf irregular galaxies (e.g., van Zee et al. 19 9 8, van Zee fc 
Haynes 2006). For example, llzotov fc ThuanI (|l999t ) found 
mean values of -0.72 ±0.06, -1.56 ±0.06, and -2.26 ±0.09 
for the log(Ne/0), log(S/0), and log(Ar/0) ratios, respec- 
tively. 

Despite their higher uncertainties, we notice that the 
Ne/O ratio seems to increase slightly with in creasing oxy- 
gen a, bundance. This effect was reported by llzotov et al.l 
(|2006l ). who interpreted it as due to a moderate deple- 
tio n of oxygen onto g rains in the most metal-rich galax- 
ies. l^rma^eOD (|2003D observed an underabundance of sul- 
phur in relatively high-metallicity starburst galaxies. These 
authors interpreted this effect as a consequence of the de- 
pletion of sulphur onto dust grains. However, we do not 
see any underabundance of sulphur; quite the opposite, the 
sulphur abundance is well correlated with both neon and 
argon abundances. 

On the other hand, it seems that the Fe/O ratio slightly 
decreases with increasing o xygen abundance. This effect 
was previously observed by llzotov et al.l (j2006D , who sug- 
gested that it may be a consequence of the iron depletion 
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Table 4. Results of the comparison between the oxygen abundance given by several empirical calibrations and the oxygen 
abundance derived here following the direct (Tc) method. 



Parameter 






R23 








N2 




03iV2 


Calibration" 


POl 


PT05 


N06 M91 


KD02 


KK04 


D02 


PP04 


N06 


PP04 N06 


Average 
Notes'* 


0.07 
0.05 

B/A (1) 


0.08 
0.07 

B/A 


0.14 0.15 
0.12 0.11 
(2) S.H. 


0.28 
0.18 
S.H. 


0.27 
0.13 
S.H. 


0.14 
0.10 
S.H. (3) 


0.12 
0.10 

B/A (4) 


0.18 
0.14 
S.L.(5) 


0.12 0.21 
0.10 0.16 
B/A (6) S.L. 



The names of the calibrations are the same as in Table \K?2\ 
^ Average value (in absolute values) of the difference between the abundance given by empirical calibrations and that obtained 
using the direct method. The names of the calibrations are the same as in Table IA.2I 

'^ Dispersion (in absolute values) of the difference between the abundance given by empirical calibrations and that obtained using 
the direct method. 

d We indicate if the empirical calibration gives results both below and above the direct value (B/A), if they are systematically 
higher than the direct value (S.H.) or if they are systematically lower than the direct value (S.L.). Some additional notes are: 

(1) Higher deviation in the low branch. 

(2) This calibration provides lower oxygen abundances in low-metallicity regions and higher oxygen abundances in 
high-metallicity regions. 

(3) Systematically higher only for 12+log(0/H)>8.2. 

(4) Higher deviation for 12+log(O/H)<8.0. Considering 12+log(O/H)>8.0, we get average=0.08 and cr=0.06. 

(5) Higher deviation at lower oxygen abundances. 

(6) Higher deviation for 12+log(O/H)<8.0. Assuming 12+log(O/H)>8.0, average=0.09 and (t=0.06. 



in dust grains, which is more important for galaxies with 
higher metallicities. 

Finally we remark how important is to use high-quality 
data and a proper estimate of all the physical parameters 
(including Tc via the direct method) in an homogeneous 
sample of objects to get reasonable conclusions about these 
topics (i.e., this work; Hagele et al. 2008). 

4.4. Comparison with empirical calibrations 

We used the data of the 31 regions for which we have a 
direct estimate of To and, hence, a direct estimate of the 
oxygen abundance, to check the reliability of several em- 
pirical calibrations. A recent review of 10 metallicity cali- 
brations, incl uding theoretical and em pirical methods, was 
presented by iKewlev &: EllisonI ()2008D . Appendix A gives 
an overview of the most common empirical calibrations and 
defines the typical parameters that are used to estimate 
the oxygen abundance following these relations. These pa- 
rameters are ratios between bright emission lines, the most 
commonly used are i?23, y, ^2, and O3N2 (see defi- 
nitions in Appendix A) . Table lA.ll lists the values of all 
these parameters derived for each region with a direct es- 
timate of the oxygen abundance (see Paper II for details) . 
Table lA.ll also includes the value derived for the q parame- 
ter (in units of c m s~^) obtained from the optimal calibra- 
tion provided by iKewlev fc Dopital ()2002[ 1. The results for 
the oxygen abundances derived for each object and empiri- 
cal calibration are listed in Table ET2I This table also indi- 
cates the branch (high or low metallicity) considered in each 
region when using the R23 parameter (see Appendix A), al- 
though, as is clearly specified in the table, for some objects 
with 8.00<12-|-log(O/II)<8.3 we assumed the average value 
found for the lower and upper branches. 

Looking at the data compiled in Table IA.2I the huge 
range of oxygen abundance found for th e same object us- 
ing di fferent calibrations is evident. As IKewlev &: EllisonI 
(|2008ll concluded, it is critical to use the same metal- 
licity calibration when comparing properties from dif- 



ferent data sets or investigate luminosity-metallicity or 
mass-metallicity relations. Furthermore, abundances de- 
rived with such strong-line methods may be significantly 
biased if the objects under study have different structural 
properties (hardness of the ionizing radiation field, mor- 
phology of the nebulae) than those used to calibrate the 
methods (!Stasihska"2009). 

Figures [T5] and plots the ten most common calibra- 
tions and their comparison with the oxygen abundance ob- 
tained using the direct method. We performed a simple 
statistic analysis of the results to quantify the goodness of 
these empirical calibrations. Table H] compiles the average 
value and the dispersion (in absolute values) of the differ- 
ence between the abundance given by empirical calibration 
and that obtained using the direct method. We check that 
the empirical cali bration that prov ides the best results is 
that proposed by iPilvuginI ()2001al lbl) , which gives oxygen 
abundances very close to the direct values (the differences 
are lower than 0.1 dex in the majority of the objects), and 
furthermore it possesses a low dispersion. We note however 
that the largest divergences found using this calibration 
are in the low-metalli city regime. The update o f this cal- 
ibration presented by iPilvugin fc Thueinl ()2005[ ) seems to 
partially solve this problem, the abundances provided by 
this calibration also agree very well with those derived fol- 
lowing the direct method . We therefore conclude that the 
iPilvugin fc ThuanI (|2005[ ) calibration is nowadays the best 
suitable method to derive the oxygen abundance of star- 
forming galaxies when auroral lines are not observed. 

On the other hand, the resul ts given b y th e 
empirical calibrations pro vided by iMcGaugh' (T99l|), 
Kcwlcy & Dopita (200^ and lKobulnick v & Kcwlcy (2004]), 
that are based on photoionization models, are system- 
atically higher than the values derived from the direct 
method. This effect is even more marked in the last two 
calibrations, which usually are between 0.2 and 0.3 dex 
higher than the expected values. These empirical calibra- 
tio ns also have a h ig;her dispersion than th a t estim ated 
for IPilvuginI (|2001aM or IPilvugin fc ThuanI (|2005ll cali- 
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Fig. 15. Comparison between the oxygen abundances derived using the direct method (Tc, always p lotted in the x — axis) 
with tho s e e stimated usin g six different em p irical calibrat ions that consider the R23 p aramete r: M91: [ McGaugh ( 1991 ); POl: 
Pilvugin (2001a b); KD02: Kewlcv fc Dopital ((20021); KK04: iKobulnickv fc KewlevI (|2C)04"l: PT05: Pilvu gin fc ThuanI (|2005i 'l: N06: 
Nagao et al.. (j2006, ) . The bottom panel of each diagram indicates the difference between empirical and direct data. 



brations. lYin et a 1 (|2007D also found high discrepancies 
when comparing the theoretical me t allicit ies using the 
theoretical models of iTremonti et ahl (l2004l ) with the Te- 
based metallicites obtained from iPilvuginI (I2001allbh and 
iPilvugin fc fhuanl (pOOSi) . 

One of the possible explanations for the different metal- 
licities obtained between the direct method and those de- 
rived from the empirical calibrations based on photoion- 
ization models are temperature fluctuations in the ion- 
ized gas. Temperature gradients or fluctuations indeed 



cause the true metallicities based on the To-method to be 
underestimated (i.e. Peimbert 1967; Stasinska 2002,2005; 
Peimbert et al. 2007). Temperature fluctuations can also 
explain our results for NGC 5253 ("Lopcz-Sa nchez et al.l 
,2002): the ionic abundances of 0++/H+ and C++/H+ 
derived from recombination lines are systematically 0.2 
- 0.3 dex higher than those determined from the direct 
method -based on the intensity ratios of collisionally ex- 
cited lines. T his abundance discrepancy has been also found 
in Galactic (jGarci'a-Roias fc EstebanI l2007t ) and other ex- 
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Fig. 16. Comparison between the oxygen abundanc es derived using the direct method (Te, in the x — axis) with those computed 
using the calibrations using the N2 parameter -002: i Denicolo. Terlevich fc T erlevichI (1200^ 1: PP04b: [Pettini fc Paeel (■2004)1 con- 
sidcring their cubic fit to N2; NOG: iNagao et al] (|2006l l- and the O.qN^ parameter following the cubic fit provided bv lNagao et alj 
(|2006i '). The bottom panel of each diagram indicates the difference between empirical and direct data. 



tragalactic (jEsteban et al.|[2009f) Hii regions and interest- 
ingly this discrepancy is in all cases of the same order as 
the differences between abundances derived from the direct 
methods and empirical calibrations based on photoioniza- 
tion models. 

The conclusion that temperature fluctuations do exist 
in the ionized gas of starburst galaxies is very important 
for the analysis of the chemical evolution of galaxies and 
the Universe. Indeed, if that is correct, the majority of 
the abundance determinations in extragalactic objects fol- 
lowing the direct method, including those provided in this 
work, have been underestimated by at least 0.2 to 0.3 dex. 
Deeper observations of a large sample of star- forming galax- 
ies -that allow us to dete ct the faint recombin ation lines, 
such as those provided bv 'Esteba n et al.l (|2009^" and more 
theoretical work -including a better understanding of the 
photoionization rnodels. such as the analysis provided by 
iKewlev fc Ellison! (|2008f )- are needed to confirm this puz- 
zling result. 

On the other hand, we ha ve checked the val idity of 
the recent relation provided by iNaeao et al.l ([200^ , which 
merely considers a cubic fit between the R23 parameter 
and the oxygen abundance. This calibration was obtained 
combining data from several large galaxy samples, the ma- 
jority from the SDSS, which includes all kinds of star- 
forming objects. As it is clearly seen in Table IA.2I and 



in Fig. \Wi the iNagao et al.l (|2006[ ) relation is not suit- 
able to derive a proper estimate of the oxygen abundance 
for the majority of the objects in our galaxy sample. In 
general, this calibration provides lower oxygen abundances 
in low-metallicity regions and higher oxygen abundances 
in high-metallicity regions. Objects located in the metal- 
licity range 8.00<12+log(O/H)<8.15 have systematically 
12+log(O/II)N06 "-^8.07 because we have to use an average 
value between the low and the high branches. Furthermore, 
ma ny of the region s do not have a formal solution to the 
Na eao et al.l (|2006D equation, such as the maximum value 
for i?23 is 8.39 at 12-f log(O/H)=8.07. We consider that the 
use of a n ionization parameter -P as introduced bylPilvuginI 
(|2001al lbh or q as followed by iKewlev fc Dopital (|2002D - is 
fundamental to obtain a real estimate of the oxygen abun- 
dance in star-forming galaxies, especially in objects showing 
strong starbursts. In the sam e sense, the direct m ethod and 
not the formulae provided bv llzotov et al.l (|2006t ) (which as- 
sumes a low-density approximation in order not to have to 
solve the statistical equilibrium equations of the O"*"^ ion) 
provides a good approximation to the actual oxygen abun- 
dance when the auroral line [O ill] A4653 is observed. 

Empirical cali brations considering a linear fit to 
the iV9 ratio dPenicolo. Terlevich fc TerlevichI 120021: 
iPettini fc Page]|l2004 give results that are systematically 
'--^O.IS dex higher that the oxygen abundances derived 
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from the direct method. The difference is higher at 
higher metahicities. We do not consider that this trend 
is a conseque n ce of com paring different objec ts: both 
iDenicolo etall (|2002D and iPettiiii fc Paeel (|2004D cahbra- 
tions are obtained using a sample of star-forming galaxies 
similar to those a nalysed in this work, many of which 
are WR galaxies. 'Denicolo et al.l (j2002f) compared the 
N2 ratio with the ionization parameter together with the 
results of photoionization models and concluded that most 
of the observed trend of N2 with the oxygen abundance 
is caused by metallicity chan g es. T he cubic fit to N2 
performed by Petti ni &: Pagell (|2004l ) better reproduces 
the oxygen abundance, especially in the intermediate- and 
high-metallicity regime (12-|-log(O/H)>8.0), where it has 
an average e rror of ^0.08 dex. H owever, the cubic fit to N2 
provided by iNaeao et al.l ((2003) gives systematically lower 
values for the oxygen abundance than those derived using 
the direct method, having an average error of ~0.18 dex. 

The empirical calibration between the oxygen 
abundance and the O3N2 parameter proposed by 
IPettini fc Pagell (|2004[ ) gives acceptable results for objects 
with 12-l-log(O/H)>8.0, with the average error ^0.1 
dex. H owever, the new relation provided by iNagao et al.l 
()2006[) involving the O3N2 parameter gives systematically 
lower values for the oxygen abundance. As we commented 
before, we consider that the Nagao et al. (2006i) calibra- 
tions are not suitable for studying galaxies with strong 
star-formation bursts. Their procedures must be taken 
with caution, galaxies with different ionization parameters, 
different chemical evolution histories, and different star 
formation histories should have different relations between 
the bright emission lines and the oxygen abundance. This 
issue is even more important when estimating the metal- 
licities of intermediate- and high-redshift galaxies, because 
the majority of their properties are highly unknown. 



5. Metallicity-luminosity relations 

The metallicity of normal disc galaxies is strongly cor- 
related with the galaxy mass. The first mass-metallicity 
relati on was found f o r irre g ular and blue compac t galax- 
ies (jLeaueux et al.l 119791 : iKinman fc DavidsonI [l98l . 
However, luminosity is often used instead of mass because 
obtaining reliable mass estimates is difficult. Rub in et al.l 
(|1984D provided the first evidence that metallicity is 
correlated with luminosity in disk galaxies. Further 
work using larger sa mples of nearby disk galaxies con- 
firme d 
19851: 



- this result (Bothu n et al.l 11984 IWvse fc Sill^ 
Skillman et all 1989; V ila-Costas fc Edmunds! Il992t 



Zaritskv. Kennicutt fc Huchrsj 11994 iGarnettI (20021) . 



Despite the huge observational effort, the origin of the 
luminosity-metallicity is still not well understood. The 
two basic ideas are (i) it represents an evolutionary 
sequence - more luminous galaxie s have processed a 
large r fraction of their raw materials ([McGaugh fc de Blo5 
Il997t iBell fc de JomI 120001: iBoselh et al.ll200lD - or (ii) it 
is related to a mass-retention sequence -more massive 
galaxies retain a larger fraction of the i r processed materia l 
(|Garnettl [200I iTremonti etHI 12004 ISalzer et al.ll2005D . 
Furthermore, other factors may play a key role in the 
variation of the metal content of a galaxy, remarking 
the quick metal enrichment that strong star-formation 
events in dwarf galaxies, such as BCDGs, may experience. 
In these objects, the freshly processed material may be 



expelled into the intergalactic medium via galactic winds 
or be mixed with the reservoirs of non-synthesized gas, in 
both cases decreasing the global metallicity of the galaxy. 

In addition, luminosity-metallicity relations are very 
useful to discern between pre-existin g dwarf galaxies and 
tidal dwarf gala xy (TDG) candidates (jPuc fc Mirabelll 19981: 
iDuc et al]l2000[ ) because these objects should have a metal- 
licit y similar to that observed in their par ent spiral galax- 
ies (jWeilbacher. Due fc AlvenslebenI l2003h and not a low- 
metallicity as it is found in dwarf objects. 

We studied the metallicity-luminosity relation using the 
data provided by our analysis. Figure [T7| plots the oxygen 
abundance vs absolute magnitude in B and NIR filters and 
including some relationships found by previous studies. In 
this figure we distinguish between galaxies (blue points) 
and TDGs candidates (red triangles) found in HCG 31 and 
Mkn 1087 groups. We also distinguish the dwarf object 
surrounding POX 4 (labeled POX 4 comp in Figures and 
Tables) because it may be another TDG, and the galaxy 
SBS 1054-1-365 because, as we will see later, its position 
in the metallicity-luminosity diagrams is quite unusual. We 
estimated the NIR absolute magnitudes for SBS 0948-1-532 
and SBS 121H-540 assuming ^ - J ~ 0.8, J - iJ - 0.3 and 
H — Ks ^0.15, that are the average values found in objects 
with properties similar to these two galaxies (see Paper I). 

The Mr-O/H diagr am includes the relation given by 
iRicher fc McCall (fT995l) for dwarf and irregular galax- 
ies (Mb > —18) extr a polat ed to high luminosities and 
the Ivan Zee fc Havne^ (|2006t ) relation found for isolated 
dwarf irregular galaxies. Both relations have a similar slope 
(—0.147 and —0.149, respectively) and intercept (5.67 and 
5.65, respectively). Our observational data have a rather 
high dispersion, but the tendency of increasing oxygen 
abundance with increasing absolute _B-magnitude is clear. 
Most of the objects fainter than AIb——18 are located above 
those relations, but many of them are TDG candidates. On 
the other hand, a substantial fraction of the brighter objects 
tend to be clearly below the metallicity-luminosity relations 
obtained by previous authors. The best linear fit to our 
data excluding the TDG candidates -it is well-known that 
they should not follow the metallicity-luminosity relation 
(|Duc fc Mirabellll998D - provides a slope of -0.099±0.019 
and an intercept of 6.37±0.37. The slope we derive for 
our galaxy sampl e is sh allow er than those provi d ed by the 
IRicher fc McGalll (I1995D and Ivan Zee fc HavnesI ()2006D re- 
lations. However, the lTremonti et al.l (|2004 ) Mr-O/H rela- 
tion for all kinds of galaxies using SDSS data (plotted with 
a blue dotted line in Fig. [T7|) show the steepest slope of all 
relations, which has a valu e of —0.185. Th a t disa gree with 
the conclusions reached by ITremonti et al.l (|2004[) , who ex- 
plained the flattening of the M—Z relation at higher masses 
because efficient galactic winds are able to remove metals 
from low-mass galaxies^ 

Some authors (ICampo s-Aguilar Moles fc Masegosal 
I1993I : iPefia fc AvaTal 11993: 'M cGaughlll994D have alreadv 
questioned the validity of the luminosity-metallicity 
relation in starbur sting galaxies. As we explained in 
previo us papers (jLopez^anchez. Esteban fc Rodriguez! 
l2004al lb[). we consider that the emission of the dominant 
young stellar population in these galaxies is increasing 
their B-luminosity, and hence the use of the standard 
metallicity-luminosity relation is not appropriate for 
starburst-dominated galaxies. Indeed, the increment of 
the _B-luminosity is moving all star-forming objects away 
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Fig. 17. Metallicity-luminosity diagrams for the objects analysed in this work. The oxygen abundance is expressed in units of 
12+log(0/H); the luminosity is expressed as the absolute magnitude in B, J, H and Ks fil ters. Red triangles represent the TDG 
candid ates found in HCG 31 (|L6pez-Sanchez. Esteban fc Rodriguezll2004ar ) and Mkn 1087 l|L6pez-Sanchez. Esteban fc Rodn'gued 
I2004bl ). A green diamond corresponds to the dwarf object surrounding POX 4, while a yellow star indicates the galaxy 
SBS 1054+365. Lin ear fits to our data are sh own with a continuous red line in the four panels. The Ms-O/H diagram includes the 
relation derived bv lRicher fc McCalll ( 1995f) fco ntinuous dark yellow line) extrapolated to high luminosities (dashed-dottcd dark 
yellow line) , the relation provided by Ivan Zee fc Haynos ( 2006) (dashed green line) , and the relation derived by Trcmonti ct aLj 
l|2004[) using SDSS data (dotted blue l ine). The Mxg-O/H diagram includes the two relations found by Salzer et al. (2003)^ 
one consi dering the [Edmunds fc Pagell ([19841 ) calibration (EP, dashed green line) and another assuming the relationships pro - 
vided by iKennicutt et al.l (|2003l l fKPG. dashed-dotted purple line) , and the relation found by iMendez de Oliveira et alT (|2006l l 
(dashed-dotted yellow line). 



-towards more negative magnitudes- from the usual 
relations valid for non-starbursting galaxies, and even 
producing -incidentally- the TDG candidates to agree 

with the relations. 

As proposed by 'Hidalgo-G amez fc OlofssonI ()1998l ) and 
reviewed by Salzer ct al. (2005|), perhaps NIR magnitudes 
are more suitable than the optical i?-magnitude to built 
metallicity-luminosity diagrams. Indeed, NIR magnitudes 
are less affected by extinction and more directly related to 
the stellar mass of the galaxy than the optical luminosi- 
ties. Furthermore, the effect of variations in star-formation 
histories and stellar populations is less pronounced in the 
NIR than in the optical. We analysed the behaviour of the 
oxygen abundance with the J, H and Ks absolute magni- 
tudes (see Fig. [T7)) . As we expected, the oxygen abundance 
increases with the luminosity. The slopes, — 0.102±0.017, 
-0.099±0.016, and -0.098±0.015 for Mj, Mh, and Mks, 
respectively, and intercepts, 6.21±0.34, 6.23±0.33, and 



6.22±0.31, of the lineal fits are remarkably similar to the 
fit parameters of our Ms-O/H relation. However, these fits 
have a better correlation coefficient (0.773, 0.779 and 0.790) 
and dispersion (0.18, 0.17 and 0.17) than thos e derived for 
the M b-O/H relation (r =0.719 and g=0.20). ISalzer et al.l 
(|2005[ ) found a notable decreasing of the rms scatter of the 
metallicity-luminosity relation between the blue and the 
NIR. 

The Mjfg-O/H di agram shown in Fig . [T71 includes the 
relations obtained by ISalzer et al.l (j2005|) -which consider 
two different empirical calib rations of their oxyg en abun- 
dance, one following the Edm unds fc Pagel (11984) method 
and an other considering t helKennicutt. Bresolin fc GarnettI 
(|2003[ ) relationships- and lMendez de Oliveira et al. (2006^ 
These relationships have a very different slope (—0.212, 
—0.195, and —0.14, respectively) and intercept (3.919, 
4.029, and 5.55, respectively) than that derived in our anal- 
ysis of the Mifg-O/H relation. However, for objects with 
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Fig. 18. Nitrogen abundance vs. luminosity for the galaxies analysed in this work. The nitrogen abundance is expressed 
in units of 12-|-log(N/H); the luminosity is expressed as the absolute magnitude in B (left) and Ks (r ight) filters. Red 
triangles represent the TDG candidates fo und in HCG 31 ()L6pez-Sanchez. Esteban fc Rodriguez! I2004al ) and Mkn 1087 
l|L6pez-Sanchez. Esteban fc Rodriguez! l2004bl ). A green diamond co rresponds to the dwarf o bject surrounding POX 4, while a 
yellow star indicates the galaxy SBS 1054+365. The relation derived [van Zee fc Havne3 l|2006l ) is plotted in the Ms-N/H relation 
with a green dashed line. The fits to our data are shown with a continuous red line. 



Mks < — 19, the relations provided bv lSalzer et al.l (|2005[ ) 
agree fairly well with the position of our data points. These 
authors found that the slopes of the metallicity-luminosity 
relations change systematically from the shortest to the 
longest wavelengths, in that the relation is steepest in the 
blue a nd more shallow in the NIR. In fact, Salze r et al.l 
(|2005D also noted that their derived A/s-O/H relation has 
a much steeper slope than those found by previously rela- 
tions (i.e., Skillman et al. 1989, Richer & McCah 1995). 

The disagreement between our data and previous 
metallicity-luminosity relations may arise because our sam- 
ple was chosen considering WR galaxies, which host strong 
starbursts, while other samples consider any star-forming 
galaxies. It may also occur because dwarf galaxies and 
normal galaxies have different relations (i.e. Tremonti et 
al. 2004). Another fact that can contribute to the afore- 
mentioned disagreement is that the oxygen abundances 
derived in previous analyses have been usually computed 
via photoionization-based empirical calibrations. As we dis- 
cussed before, they seem to overestimate the actual O/H 
ratio by about 0.2 dex. 

In Fig [17] we can see that the separation between the po- 
sitions of TDG candidates and dwarf galaxies is somewhat 
more evident in the metallicity-luminosity diagrams involv- 
ing NIR magnitudes. POX 4 comp agrees with the relations, 
and hence further analysis of this object is needed to under- 
stand its real nature. On the other hand, SBS 1054-f 365 has 
a higher oxygen abundance for its optical and NIR absolute 
luminosities. We do not consider the O/H to be overesti- 
mated because of the high quality of our spectra (Paper II) 
and because previous studies derived the same metallicity 
(i.e., Izotov & Thuan 1998). We are also quite confident 
of the correct determination of its optical and NIR magni- 
tudes because of the detailed analysis we already presented 
in Paper I. Hence we consider that this nearby (8 Mpc) 
BCDG may have experienced a recent a strong pollution of 
heavy elements which probably has not had enough time to 
disperse and mix with the surrounding ISM. In Paper V we 
will see that this galaxy also possesses a very high Mdyn/L b 



ratio, suggesting that the neutral gas is quite disturbed. 
Further analysis of this intriguing object, including H i ra- 
dio data and optical spectroscopy, would be necessary to 
ascertain its true nature. 

Finally, another aspect that should also be addressed is 
the possibility that some of the sample objects may be a 
recent merger of two galaxies. In this case, the integrated 
luminosity would be higher than that expected for a single 
galaxy, and hence its position would be also far from the 
metallicity-luminosity relations. That may be happening in 
UM 420 -see Sect. 3.9 in Paper I, Sect. 3.9.3 in Paper II and 
iJames et al] (|2009D - and Tol 1457-262 (see Sect. 3.18.3 in 
Paper II), which show a considerable high luminosity given 
their low oxygen abundances. Hence it is more appropriate 
to compare the stellar mass (and not the B or the NIR 
luminosities) with the metallicity to reach any conclusions 
about this topic, as we will explain in Paper V. 

The nitrogen abundance also has a strong correlation 
with the luminosity. Figure [T8| plots the nitrogen abun- 
dance, expressed in units of 12-|-log(N/H), as a function 
of the absolute B and Ks magnitudes. As in Fig. [T71 we 
distinguish the TDG candidates and the intriguing objects 
POX 4 comp and SBS 1054+365. From Fig.[lH|it is also evi- 
dent that the nitrogen abundance increases with increasing 
luminosity. We performed a linear fit to our data (with- 
out considering the problematic objects already discussed) , 
which is plotted with a continuous red line. In the Mb- 
N/H diagram this fit yields a slope of —0.181 ± 0.036 and 
an intercept of 3. 65±0.68. This resul t is vir tually identical 
to that found by Ivan Zee fc Havnei (|2006[ ) , their relation 
is shown with a green dashed line in the left diagram of 
Fig.[TH| The scatter of the Afs~N/H relation is higher than 
that found in the Mb-O/H relation, actually, we may dis- 
tinguish between two kinds of objects in this figure, because 
the more evolved galaxies (such as Mkn 1199, Mkn 1087, 
Tol 9, ESO 566-8) and the non-starbursting systems 
(Mkn 1199 NE, IRAS 08339+6517 c, ESO 566-7) have a 
higher N/H value at a given absolute i?-luminosity than the 
very starbursting (i.e. IRAS 08339+6517, HCG 31 AG) and 
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Fig. 19. Met allicity- colour diagrams for the galaxies analysed in this work. A linear fit is shown with a continuous red line. 



blue compact dwarf (i.e. UM 420, POX 4, SBS 0926+606 A) 
galaxies. TDG candidates do not lie far from this rela- 
tion. These differences are probably a consequence of the 
different star-formation histories of the galaxies, as we 
also explained the scatter in the observed N/0 ratio at 
12-hlog(0/H)>7.9 (see Fig.E]). 

As it happened in the previous diagrams involving 
the oxygen abundance, the luminosity-metallicity diagram 
using NIR data has a smaller scatter (Fig. [181 right). 
The fit to our data in the Mi^g-N/H diagram, neglecting 
TDG candidates and problematic objects, gives a slope of 
-0.205 ± 0.030 and an intercept of 2.86±0.61. The correla- 
tion coefficient is r=0.813 and the dispersion is ct=0.33. 
In comparison, the Ms-N/H diagram has r^O.TO? and 
(7=0.40. We note that again TDG candidates and the 
galaxy SBS 1054-f 365 lie far from this fit, which is plotted 
in Fig. [15] (right) with a continuous red line. POX 4 comp 
agrees quite nicely with the relation. In both diagrams, the 
position of Tol 1457-262 is far from the relation, having 
a considerable high luminosity for its N/H ratio. As we 
already discussed, we consider that this object may cor- 
respond to a recent merger of two independent galaxies. 
Hence we conclude that luminosity-metallicity relations us- 
ing NIR band magnitudes are very suitable to understand 
galaxy properties and evolution, providing tighter correla- 
tions than those relations involving optical magnitudes. 



6. Metallicity-colour relations 

Wc analysed the relationship between the optical and 
NIR colours observed in each galaxy and its metallicity. 
Figure [T9l plots the results for all colours (from Table[T]) vs. 
the oxygen abundance and a fit to the data. Although the 
dispersion is large, the general trend is that galaxies with 



redder colours have a higher oxygen abundance in their 
ionized gas. This suggests again the relative importance of 
the evolved stellar populations existing in the galaxies, and 
agrees both with models of galaxy evolution (i.e., Leitherer 
et al. 1999; Bruzual & Charlotte 2003) and previous ob- 
servational analysis (i.e. Jansen et al. 2000; Lilli, Carollo 
& Stockton 2003). Neglecting other effects (such as ex- 
tinction), the birth of new generations of stars within a 
galaxy will continuously increase the number of intermedi- 
ate and low-mass stars, which are the typical stellar pop- 
ulation that constitutes the low-luminosity component in 
starburst galaxies. Hence it should be expected that as the 
system evolves -increasing its metallicity- its optical and 
NIR colours are more dominated by the evolved population, 
and therefore they are redder. Another interesting comment 
regarding these correlations is to remark how important it 
is to have an estimate of the metallicity of the galaxies when 
comparing the optical / NIR colours with the predictions of 
theoretical evolutionary models (see Paper I). 

Besides the N/0 ratio has a large dispersion in our sam- 
ple objects, we checked whether a correlation exists between 
this ratio and the B — V colour. Figure [201 plots both quan- 
tities considering the B — V colour of the complete galaxy 
(left panel) and that derived for the old underlying popula- 
tion (right panel). Both diagrams include a fit to our data 
(continuous red line) . Figure [20] shows the large dispersion 
of the data points, but there seems to be a suggestion of a 
tendency of increasing N/0 with increasing B — V colour. 
This tendency does not become clearer when considering 
the data of the underlying stellar component. However, the 
rel ation found using our dat a disagrees with that obtained 
by Ivan Zee fc Havnei (j2006[ ) analysing a sample of dwarf 
irregular galaxies. We note that they referred to the B ~V 
colour of the underlying population, but actually it is just 
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— Linear fit to our data: log (N/0) = -1 .32+0.52*(B-V) 
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- Linear fit to our data: log (N/0) = -1 .63+1 .23*(B-V)^j^ 

- VZH06: log (N/0) = -1 .96+1 .26*(B-V) 



Fig. 20. N/0 ratio vs. the reddening-corrected B — V colour considering all the galaxy (left) or just the underlying stellar 
population (right). A linear fit to our data is shown with a continuous red line. The relationship found bv ivan Zee &z. Havne^ 
l(2006f) in their analysis of dwarf irregular galaxies is plotted with a dashed green line. 
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Linear fit to ttie data: log (N/0)= -1 .42+0.55*(J-H) 
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Fig. 21. N/0 ratio vs. the reddening-corrected V ~ R (left) and J — H (right) colours. A linear fit to our data is shown with a 
continuous red line. 



the integrated colour of the galaxy. The significant differ- 
ence found between both samples may be a consequence 
of two reasons: (i) we are comparing different objects -our 
sample is co mposed by starburst i ng gal axies, many of them 
BCDGs, but Ivan Zee fc HavnesI ()2006[) analysed quiescent 
dwarf irregular galaxies with star-formation histories that 
are correctly modelled by simple changes in the SFR- and 
(ii) the N/O rat ios tabulated by these authors, who used 
lMcGaugh[(jl99lt) photoionization models in many cases, are 
not well estimated -as we saw before, these models over- 
estimate the oxygen abundance, and hence underestimate 
the N/O ratio, as it is clearly seen in the 0/H-N/O di- 
agram (Fig. [TTl right) when comparing the Ivan Zee et al.l 
(|1998() sample with the results provided by other authors-. 
We consider that the first reason may be the most rea- 
sonable explanation: objects with different star-formation 
histories have different relationships between the N/O ratio 
and other global characteristics such as optical colours. 

However, as we commented before, using the B- 
magnitude may be not ideal to compare the colours and the 
chemical properties of these galaxies. Hence we also com- 
pared the N/O ratio with other optical and NIR colours, 



two of them {V — R and J — H) are shown in Fig. j^H As we 
already explained when analysing the N/H-L relationship 
(Fig. [TS|) . there are two kinds of objects with similar op- 
tical luminosities but with different nitrogen abundances. 
This is also reflected in the N/O ratio (Fig. ^KA left), where 
normal and dwarf non-starbursting galaxies show a higher 
N/O ratio than starbursts (which have a lower N/O ratio in 
Fig. HOI left). This segregation is not so evident when com- 
paring the N/O ratio with the NIR luminosities (Fig. 
right), perhaps because the underlying component makes a 
higher contribution to the NIR magnitudes, but more data 
are needed to quantify this effect. 

Finally, we also note that we do not see any clear rela- 
tionship between the underlying population colours and the 
N/O ratio: for example, the V—R colour of the low luminos- 
ity component is 0.31 ±0.05 for all ranges of N/O. It was not 
possible to analyse the NIR colours of the low-luminosity 
component in the majority of the cases (see Table 5 in 
Paper I), but for the few available data it also seems that 
they are metallicity-independent, showing J — H ^ 0.5±0.1 
and H -Ks ^ 0.2. 
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7. Conclusions 

We compiled and analysed globally the optical/ NIR colours 
and the physical and chemical properties of the ionized 
gas for a sample of 20 Wolf-Rayet galaxies. The individual 
analysis of the photometry of each galaxy was presented 
in Paper I, while the individual spectroscopic analysis was 
discussed in Paper II. The metallicity of these galaxies lies 
between 7.58 and 8.75 -in units of 12-|-log(0/H)-. The most 
important conclusions found in this study are 

1. The colours estimated for our galaxy sample, which were 
corrected for both extinction and nebular emission us- 
ing our spectroscopic data, agree quite well with the 
predictions given by evolutionary synthesis models, es- 
pecially in compact and dwarf objects. Small discrep- 
ancies are explained because of the existence of several 
stellar populations within each galaxy and differences in 
their star-formation history. All galaxies show evidence 
of an old stellar population underlying the starburst, 
with ages older than 500 Myr. 

2. We checked that all objects can be classified as pure 
star-forming galaxies. In total, we compiled 41 indepen- 
dent star-forming regions, of which 31 have a direct es- 
timate of the electron temperature of the ionized gas, 
and hence their element abundances were derived using 
the direct method. We found that younger bursts have 
larger ionization budgets and are therefore capable to 
heat the ionized gas to higher electron temperatures. 
Both c(II/3) and Wabs increase with increasing metallic- 
ity, as predicted by galaxy evolution models. 

3. We compiled the oxygen abundance and N/0, S/0, 
Ne/0, Ar/0, and Fe/0 ratios in our WR galaxy sample. 
They generally agree well with previous observations. 
The N/0 ratio is found to be rather constant for objects 
with 12-|-log(0/H)<7.6, has an important dispersion 
in galaxies with 7.6<12-|-log(0/II)<8.3, and increases 
with the metallicity in objects with 12-|-log(0/H)>8.3. 
This behaviour is explained assuming the very different 
star-formation histories that each individual system has 
experienced. 

4. We detected a high N/0 ratio in objects showing clear 
WR features (HCG 31 AC, UM 420, IRAS 0828+2816, 
III Zw 107, ESO 566-8 and NGC 5253). The ejecta of 
the WR stars may be the origin of the N enrichment in 
these galaxies, but further detailed data comparing the 
chemical properties of a larger sample of both WR and 
non-WR galaxies, as well as careful analyses of galaxies 
showing a localized high N/0 ratio, are needed. 

5. The relative abundance ratios of the a-elements to oxy- 
gen are approximately constant, which is expected be- 
cause all four elements are mainly produced by mas- 
sive stars. We found indications of a moderate depletion 
of oxygen and iron onto grains in the most metal-rich 
galaxies. 

6. We compared the abundances provided by the direct 
method with those obtained using the most common 
empirical calibrations: 

- The Pilyugin- method of iPilvueinl (|2001al lbh. which 
considers the R-2?, and the P para. meters and is up- 
dated bv lPilvugin k, ThuanI ()2005[ ). is nowadays the 
best suitable empirical calibration to derive the oxy- 
gen abundance of st ar-forming ga l axies. The cubic fit 
to i?23 provided by iNagao et al.l (|2006D is not valid 
for analysing these star-forming galaxies. 



The relations between the oxygen abundance and 
the N2 or the O^N 2 parameters provided by 
iPettini fc Pagell ([2004) give acceptable results for 
objects with 12-}-log(O/H)>8.0. 
The results provided by empirical calibrations 



based on photoiqnization models (iMcGaugh 19911 



iKewlev fc Dopitall2002HKobulnickv fc Kewlevl200- 

are systematically 0.2 - 0.3 dex higher than the 
values derived from the direct method. These dif- 
ferences are of the same order as the abundance 
discrepancy found between abundances determined 
from recombination and collisionally excited lines 
of heavy-element ions. This may suggest tempera- 
ture fluctuations in the ionized g they exist in 
Galactic and other extragalactic H 11 regions. 

7. We studied the optical/ A^ZR luminosity-metallicity re- 
lations for our sample galaxies. We found that our data 
generally disagree with previous relations, perhaps be- 
cause the objects analysed here host strong starbursts, 
but maybe also because we used the direct method 
and not the empirical calibrations to derive the oxy- 
gen abundances. The L-Z relations tend to be tighter 
when using NIR luminosities, TDG candidates are in- 
deed easier detected using the Mks-Z relation. 

8. The nitrogen abundance also has a strong correla- 
tion with the luminosity, but normal and dwarf non- 
starbursting galaxies show a higher N/0 ratio than 
strong starbursting galaxies. 

9. We found that galaxies with redder colours tend to 
have a higher oxygen abundance. The N/0 ratio also 
increases with redder colours. Both results agree with 
galaxy evolution models. The colours of the underly- 
ing component seem to be metallicity-independent, but 
more data are still needed to confirm this trend. 

We finally conclude that it is fundamental to perform a de- 
tailed analysis of both the photometric (optical and NIR 
magnitudes and colours of both burst and underlying com- 
ponent) and the chemical (oxygen abundances, N/0 ratios) 
properties of these star-forming galaxies to understand the 
evolutionary stage of every system. A larger galaxy sam- 
ple following the ideas compiled in this work will comple- 
ment the results derived from huge databases (which do 
not distinguish the evolutionary state of the galaxies or 
the relative contribution of the burst /underlying popula- 
tions, and for which the majority of the properties have 
been derived automatically and, in many cases, have consid- 
ered empirical calibrations to determine the chemical abun- 
dances) when analysing the global properties and evolution 
of star- forming galaxies. 
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Appendix A: Empirical calibrations of the oxygen abundance 

When the spectrum of an extragalactic H ii region does not show the [O ill] A4363 emission line or other auroral lines 
that can be used to derive Tc, the so-called empirical calibrations are applied to get a rough estimation of its metallicity. 
Empirical calibrations are inspired partly by photo-ionization models and partly by observational trends of line strengths 
with galactocentric distance in gas-rich spirals, which are believed to be due to a radial abundance gradient with abun- 
dances decreasing outwards. In e xtragalactic objects, the usefulness of the empirical methods goes beyond the derivation 
of abundance gradients in spirals (jPilvugin. Vilchez fc Continill2004[) , as these methods find application in chemical a bun- 
dance studies of a variety of objects, including low-surface brightness galaxies (|de Narav. McGaugh Sz de Blokll2004f) and 
star-forming galaxies at intermediate and high redshift, where the advent of 8-10 m class telescopes has made it possible 
to extend observations (e.g., Teplitz et al. 2000, Pettini et al. 2001; Kobulnicky et al 2003; Lilly, Carollo & Stockton 
2003; Steidel et al. 2004; Kobulnicky & Kewley 2004; Erb et al. 2006). 

As the brightest metallic lines observed in spectra of H II regions are those involving oxygen, this element has been 
extensively used to get a suitable empirical calibration. Oxygen abundance is important as one of the fundamental 
characteristics of a galaxy: its radial distribution is combined with radial distributions of gas and star surface mass 
densities to constrain models of chemical evolution. Parameters defined in empirical calibrations evolving bright oxygen 
lines are 

/([OIII])A4959 + /([OIII])A5007 ^, 
i?3 = , (A.l) 

_ /([OII])A3727 

' ^^-2) 

R23^R3 + R2, (A.3) 

^ = 1^, (A.4) 

^23 

i?3 1 

iJensen. Strom fc StromI (jl976|) presented the first explorati on in this method considering the R^ index, which con- 
siders the [Oili] AA4959,5007 emission lines. However, were iPagel et al.l (|1979| ) who introduced the most widely used 
abundance indicator, the R23 index, which also included the bright [ Oil] A37 27 emission line. Since then, many stud- 
ies have been performed to refine the calibration of R23 (Edmunds fc Pagel 1984; McCall. Rvbski & Shields. 198^ 



Dopita fc Ev ans 1986; Torres-Peimber t et al.l 119 89: McGaugh 1991; Zaritskv, Kcnnicutt fc Huchra 1994; Pilvugin 200C 
2001allbl : lKewlev fc Dopitaii2002i:iKobulnickv fc Kewlevii2004; Pilvugin fc Thuaa200 5: Nagao. Maiolino fc Marconi„200e '" 



The most successful are the calibrations of lMcGaughl (Il99lh an d Kc wlev fc Dopita|_( 2002'). which are based on photoion- 
ization models, and the empirical relations provided by FilvuginI (|2001allbf ) and iPilvugin fc ThuanI ()2005[ ). Both kinds 
of calibrations improve the accuracy by making use of the [O ill] /[On] ratio as ionization parameter, which accounts 
for the large scatter found in the R23 v ersus oxygen abundance calibration, which is larger than observational errors 
(|Kobulnickv. Kennicutt fc Pizagnolll999D . Figure \A.1\ shows the main empirical calibrations that use the R23 parameter. 

The main problem associated with the use of R23 parameter is that it is bivaluated, i.e., a single value of R23 can be 
caused by two very different oxygen abundances. The reason of this behaviour is that the intensity of oxygen lines does 
not indefinitely increase with metallicity. Thus, there are two branches for each empirical calibration (see Fig. lA.lT) : the 
low-metallicity regime, with 12-|-log(0/H)<8.1, and the high-metallicity regime, with 12-|-log(0/H)>8.4. That means 
that a very large fraction of the star-forming regions lie in the ill-defined turning zone around 12-|-log(O/H)~8.20, where 
regions with the same R23 value have oxygen abundances that differ by almost an order of magnitude. Hence, additional 
information, such as the [Nil] /Ha or the [Oii]/[Oiii] ratios, is needed to break the degeneracy between the high and low 
branches (i.e., Kewley fc Dopita, 2002). Besides, the R23 method requires that spectrophotometric data are corrected by 
reddening, which effect is crucial because [0 11] and [O ill] lines have a considerably separation in wavelength. 

Here we list all empirical calibrations that were considered in this work, compiling the equations needed to derive the 

oxygen abundance from bright emission line ratios following every method. 

Edmund & Pagel (1 984): Although th e R23 p arameter was firstly proposed bv lPagel et al.l (|l979r ). the first empirical 
calibration was given bv lEdmunds fc Pagell ([l98^, 

12 + log(0/H)„p = 8.76 - 0.69 logi?23, 

12 + log(0/H)z„„ = 6.43 + 1.67 log i?23, (A.6) 

with the limit between the lower and the upper branches at 12-|-log(O/H)'^8.0. 

McCall, Rybski &: Shields (1985) presented an empirical calibration for oxygen abundance using the R23 param- 
eter, only valid for 12+log(0/H)>8.15. However, they did not give an analytic formulae but only listed it numerically 
(see their Table 15). The four-order polynomical fit for their values gives the following relation: 

12 + log(0/H)„p = 9.32546 - 0.360465a; -f 0.203494a;2 -1-0.278702x3 - 1.36351a;^ (A.7) 
with X = log R23 ■ 
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Fig. A.l. Empirical calibrations of oxygen abundance using the R23 parameter. Note that they are bi-valuated. The dashed zone 
indicates the region with higher uncertainties in O/H. The empirical calibrations plotted in the figure are: EP94: Edmund & Pagel 
(1984); M91: McGaugh (1991) using y^O {R2 = R3); POl: Pilyugin (2001) using P = 0.5 {R2 = -R3); (KD02-I-KKQ4): Kewley & 
Dopita (2002) using the formulation of Kobulnicky & Kewley (2004) assuming g = 7.5 x lO'^ cm s"^ N06: Nagao et al. (2006) 
using their cubic fit to J?23. 



Zaritzky, Kennicutt &: Huchra (1994) provided a simple analytic relation between oxygen abundance and -R23: 

12 + log(0/H)„p = 9.265 - 0.33a; - 0.202x2 - 0.207a;3 - 0.333a;^ (A.8) 

Their formula is an average of three previous calibrations: lEdmunds fc Pagell ()1984f ). McCall et al. (1985) and Dopita & 
Evans (1986). Following the authors, this calibration is only suitable for 12+log(O/H)>8.20, but perhaps a more realistic 
lower limit is 8.35. 

McGaugh (1991) calibrated the relationship between t he R23 ratio and g as-phase oxygen abundance using Hii 
region models derived from the photoionization co de Cloudy (| Ferland et al.l ll998 1. McGaugh's models include the effects 



of dus t and variations in ionization parameter, v. iKobulnickv et al. (1999) give analytical expressions for the McGaugh 
(|1991h cahbration based on fits to photoionization models; the middle point between iDoth branches is 12+log(0/H)^8.4: 

12 + log(0/H)„p = 7.056 + 0.767a; + 0.602a;2 - ?;(0.29 + 0.332a; - 0.331a;2), (A.9) 
12 + log(0/H)/o„ = 9.061 - 0.2a; - 0.237a;2 - 0.305a;3 0.0283a;'* 

-2/(0.0047 - 0.0221a; - 0.102a;2 - 0.0817a;3 - 0.00717a;4). (A.IO) 

Pilyugin (2000) found that the previous calibrations using the R23 parameter had a systematic error depending on 
the hardness of the ionizing radiation, suggesting that the excitation parameter, P, is a good indicator of it. In several 
papers, Pilyugin performed a detailed analysis of the observational data combined with photoioniz ation m odels to obtain 
empirical calibrations for the oxygen abundance. iPilvu gin ( 2000i) confirmed the idea of . McGaughl (|199lD that the strong 
lines of [0 11] and [O III] contain the necessary information for the determination of accurate abundances in low-metallicity 
(and may be also in high-metallicity) H 11 regions. He used new observational data to propose a linear fit involving only 
the i?23 parameter, 

12 + log(0/H)„p = 9.50 - 1.40 log R23, (A.U) 
12 + log(0/H)j„^ = 6.53 + 1.40 log i?23, (A.12) 

assuming a limit of 12 -|-log(O/H)^8.0 between th e two b r anche s. This calibration is close to that given by 
lEdmunds fc Pagel ()1984D ; it has the same slope, but iPilvuginI ()2000() is shifted towards lower abundances by around 
0.07 dex. However, this new relation is not sufficient to explain the wide spread of observational data. Thus, Pilyugin 
(2001a) give the following, more real and complex, calibration involving also the excitation parameter P: 

.0 , /^/TTX i?23 + 54.2 59.45P + 7.31P2 

12 + log(0/H)„p = 6.oi + 6.71P-H0.371P^ + 0.243P23 - ^^"'^^ 
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This is the so-called P-method, which can be used in moderately high-metallicity H ii regions (12-|-1or(0/H)^ 8.3). Pilyugin 
used two-zone models of Hii regions and assumed the re(Oii) - Te(Oiii) relation from iGarnettI (jl992[ ). For the low 
metallicity branch, Pilyugin (2001b) found that 

12 + log(0/H),„,„ = 6.35 + 1.45 log i?23 - 1.74 log F. (A.14) 

Pilyugin estimates that the precision of oxygen abundance determination with this method is around 0.1 dex. 

Pilyugin & Thuan (2005) revisited these calibrations including more spectroscopic measurements of H ii regions 
in spiral and irregular galaxies with a measured intensity of the [O ill] A4363 line and recalibrate the relation between 
the oxygen abundance and the R23 and P parameters, yielding to: 

10^1 rn/m i?23 + 106.4 + 106.8P-3.40P^ 

12 + log(0/H),.. = 17.72 + 6.60P + 6.95P^- 0.302^23 ' ^^"'^^ 
10^1 ^n/m P23 + 726.1 + 842.2P + 337.5P^ 

12 + log(0/H)„, = 85.96 + 82.76P + 43.98P2 + 1.793P23- ^^"'^^ 

Kewley & Dopita (2002) used a combination of stellar population synthesis and photoionization models to develop 
a set of ionization parameters and abundance diagnostic based only on the strong optical emission lines. Their optimal 
method uses ratios of [Nil], [On], [Oiii], [Sii], [Siii] and Balmer lines, which is the full complement of strong nebular 
lines accessible from the ground. T hey also recommend proc edures for the derivation of abundances in cases where only 
a subset of these lines is available. iKewlev &: Dopital ([2002) models start with the assumption that i?23, and many of 
the other emission-line abundance diagnostics, also depends on the ionizat ion parameter q = c x U , that has units 
ofcm s~^. They used the stellar population synthesis c odes Starburst 99 (jLeitherer et al.lll999t IVazauez fc Leithereil 
l2005[l and Pegase.2 (jPioc &: Rocca-Volmerang3ll997f) to generate the ionizing radiation field, assuming burst mod- 
els at zero age with a Salpeter IMP and lower and upper mass limits of 0.1 and 120 Mq, respectively, with metal- 
licities between 0. 05 and 3 times solar. The ionizing radiation fields were input into the photoionization and shock 
code. Mappings ( Sutherland fc Dopital 119931 ). which includes self-consistent treatment of nebular and dust physics. 
iKewlev &: Dopita ( 20021) pre viously used these models to simulate the emission- line spectra of H 11 regions and starburst 
galaxies (|Dopita et aTEoOOD . and are completely described in their s tudy. 



Kobulnicky & Ke wley (2004 ) gave a parameter ization of the Kewlev fc Dopital ()2002|) i?23 method with a form 



similar to that given by iMcGaughl ()1991[ ) calibration. iKobulnickv &: KewlevI (|2004D presented an iterative scheme to 



resolve for both the ionization parameter q and the oxygen abundance using only [Oiii], [On] and H/3 lines. The 
parameterization they give for q is 

32.81 - 1.153J/2 + [12 -I- log(0/H)] [ - 3.396 - 0.025y -I- 0.1444y2] 

~ 4.603 - 0.3119y - 0.163y2 + [12 + log(0/H)] [ - 0.48 + 0.0271y + 0.02037y2] ' ^ " 

where y = log([0 iii]/[0 11]). This equation is only valid for ionization parameters between 5x10^ and 1.5x10^ cm s^^. 
The oxygen abundance is parameterized by 

12 + log(0/H)„p = 9.72 - 0.777a; - 0.951x^ ~ 0m2x^ ~ O.Sllx^ - log(g) 

X (0.0737 - 0.0713a; - O.Mla;^ + 0.0373a;^ - 0.058a;^), (A.18) 
12 + log(0/H)/oi„ = 9.40 -t- 4.65a; - 3.17x^ - log(9)(0.272 + 0.547a; - 0.513x^), (A.19) 

being x — logP23- The first equation is valid for 12+log(0/H)>8.4, while the second for 12-|-log(0/H)<8.4. Typically, 
between two and three iterations are required to reach convergence. Following the authors, this parameterization should 
be regarded as an improved , implementation-friendly approach to be preferred over the tabulated P23 coefficients given 
by IKewlev fc Dopital (|2002D . 

Nagao, Maiolino &: Marcani (2006) did not consider any ionization parameter. They merely used data of a large 
sample of galaxies from the SDSS to derive a cubic fit to the relation between i?23 and the oxygen abundance, 

logP23 = 1-2299 - 4.1926y-h 1.0246y2 - 0.063169y^ (A.20) 

with ?;=12+log(0/H). 

Besides P23, additional parameters have been used to derive metallicities in star- forming galaxies. Without other 
emission lines, the N2 parameter, which is defined by 

/(|N11|)A6.83 ^ 

Ha 

can be used as a crude estimator of metallicity. However, we note that the [Nil] /Ha ratio is particularly sen- 
sitive to shock excitation or a hard radiation field from an AGN. The N2 parameter was firstly suggested by 
IStorchi-Bergmann. Calzetti fc KinnevI (|1994[ ). who gave a tentative calibration of the oxygen abundance using this pa - 
rameter. This ca l ibrati on h as been revisited byFvan" Zee. Salzer fc HavnesI (|1998l) ; lDenicol6. Terlevich fc TerlevichI (|2002D : 
iPettini fc Pagell (|20Q4D and lNagao et all (j2006l ). The Denicolo et al. (2002) calibration is 

12 + log(0/H) = 9.12 + O.73N2, (A.22) 
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which considerably improves the previous relations because of the inclusion of an extensible sample of nearby extragalactic 
Hii regions. The uncertainty of this method is '^0.2 dex because N2 is sensitive to ionization a nd 0/N variations, 
so strictly speaking it should be used mainly as an indicator of galaxy-wide abundances. Dcnicol d et al.l (2002) also 
compared their method with photoionization models, concluding that the observed A'2 is consistent with nitrogen being 
a c ombination of both pri mary and secondary origin. 

iPettini fc Pagell ()2004|) revisited the relation between the N2 parameter and the oxygen abundance including new 
data for the high- and low-metallicity regimen. They only considered those extragalactic H 11 regions where the oxygen 
values are determined either via the Te method or with detailed photoionization modelling. Their linear fit to their data 
is 

12 + log(0/H) = 8.90 + 0.57iV2, (A.23) 

which has both a lower slope and zero-point that the fit given bv lDenicolo et all (|2002[ ). A somewhat better relation is 
provided by a third-order polynomical fit of the form 

12 -I- log(0/H) 9.37 + 2.0327V2 + 1.26(7V2)^ + 0.32(7V2)^ (A.24) 

valid in the range —2.5 < N2 < — 0.3. lNagao et ahl ()2006() also provided a relation between N2 and the oxygen abundance, 
their cubic fit to their SDSS data yields 

log N2 = 96.641 - 39.9412/ + b.2221y^ - 0.22040y^ (A.25) 
with ?/=12+log(0/H). 

Pettini &: Pagel (2004) revived the O3N2 parameter, previously introduced bv lAlloin et al.l (|1979f) and defined by 
^ , [OIII] A5007/H/? 

^-^-"^°g [NII]A6583/Ha - ^^'^^^ 
iPettini fc Pagell (|2004l ) derived the following least-square linear fit to their data: 

12 + log(0/H) = 8.73 - O.32O37V2. (A.27) 

iNagao et al.l ()2006l ) also revisited this calibration and derived a cubic fit between the O3N2 parameter and the oxygen 
abundance, 

log 03iV2 = -232.18 + 84.423y - 9.9330?/^ + 0.37941?/^ (A.28) 
with 2/=12-flog(0/H). 

Other important empirica l calibrations that were n ot used in this study involv e the 6*23 parameter, intro- 
duced bv Ivilchez fc EstebanI (jl996[ ) and revisited by iDiaz fc Perez-Montercl (j2000t) : lOev fc Shieldi (|2000[ ) and 



iPerez-Montero fc Diad (|2005l ). In the last years, bright emissio n line ratios such as [Ar ill]/[0 ill] and [Sili]/[Oiil 



(Stasihska 2006) or [Neiii]/[Oiii] and [Oiii]/[Oii] (Nagao e t aLll2006|) ha ve been explored as indicators of the oxy 



gen abundance in H 11 regions and starburst galaxies. P eimbert et al.l (2007j suggested to use the oxygen recombination 
lines to get a more precise estimation of the oxygen abundance. Nowadays, there is still a lot of observational and the- 
oretical work to do involving empirical calibrations (see recent review by Kewley & Ellison 2008), but these methods 
shou ld be used only for objects whose H il regions have the same structural properties as those of the calibrating samples 
(fStasihska 2009 ). 
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Table A.l. List of the parameters used to compute the oxygen abundance in all regions with a direct estimation of To using 
empirical calibrations. 



Region fl2a P = Rz/R23 y =\og{R-i / R2) N2 O3N2 <?Kg02o° 



HCG 31 AC 


5.42 


0, 


.571 


0.125 


0, 


,104 


1, 


.349 


3, 


,76E+07 


HCG 31 B 


7.93 


0, 


.408 


-0.162 


0, 


.101 


1, 


,381 


4, 


,91E+07 


HCG 31 E 


7.12 


0, 


.511 


0.020 


0, 


.090 


1, 


,486 


7, 


,40E+07 


HCG 31 Fl 


8.91 


0, 


.819 


0.656 


0, 


,034 


2, 


,201 


5, 


,78E+07 


HCG 31 F2 


7.60 


0, 


.724 


0.418 


0, 


,036 


2, 


,064 


6, 


,28E+07 


HCG 31 G 


8.20 


0, 


.499 


-0.002 


0, 


,106 


1 


,462 


6, 


,96E+07 


Mkn 1199 C 


7.69 


0, 


.809 


0.627 


0, 


,131 


1, 


,555 


1, 


,55E+08 


Haro 15 A 


9.73 


0, 


.884 


0.881 


0, 


,027 


2, 


.378 


8, 


,55E+07 


Mkn 5 Al 


7.58 


0, 


.748 


0.473 


0, 


,051 


1, 


.915 


6, 


,96E+07 


Mkn 5 A2 


8.19 


0, 


.702 


0.372 


0, 


,049 


1, 


.944 


1, 


,72E+08 


IRAS 08208+2816 C 


7.77 


0, 


.793 


0.583 


0, 


,129 


1, 


,558 


8, 


,55E+07 


POX 4 


10.68 


0, 


.906 


0.986 


0, 


,015 


2, 


,697 


1, 


,05E+08 


UM 420 


6.45 


0, 


.649 


0.268 


0, 


,099 


1, 


,497 


4, 


,81E+07 


SBS 0926+606A 


7.40 


0, 


.811 


0.632 


0, 


,026 


2, 


,227 


6, 


,68E+07 


SBS 0948+532 


8.85 


0, 


.874 


0.843 


0, 


,022 


2, 


,430 


2, 


,54E+08 


SBS 1054+365 


9.33 


0, 


.893 


0.920 


0, 


,020 


2, 


,503 


9, 


,10E+07 


SBS 1211+540 


7.22 


0, 


.892 


0.918 


0, 


,008 


2, 


.788 


1, 


,16E+08 


SBS 1319+579A 


9.92 


0, 


.908 


0.996 


0, 


,014 


2, 


.671 


1, 


,05E+08 


SBS 1319+579B 


7.13 


0, 


.722 


0.415 


0, 


,046 


1, 


,922 


6, 


,15E+07 


SBS 1319+579C 


7.11 


0, 


.710 


0.389 


0, 


,052 


1, 


,860 


5, 


,91E+07 


SBS 1415+437C 


5.22 


0, 


.783 


0.558 


0, 


,015 


2, 


,301 


5, 


,91E+07 


SBS 1415+437A 


4.86 


0, 


.810 


0.629 


0, 


,012 


2, 


.370 


5, 


,44E+07 


III Zw 107 A 


7.13 


0, 


.701 


0.369 


0, 


,100 


1, 


.573 


5, 


,78E+07 


Tol 9 INT 


4.58 


0, 


.689 


0.345 


0, 


,252 





.973 


4, 


,16E+07 


Tol 9 NOT 


4.78 


0, 


.629 


0.230 


0, 


,287 





.894 


3, 


,39E+07 


Tol 1457-262A 


9.89 


0, 


.773 


0.532 


0, 


,033 


2, 


.236 


9, 


,91E+07 


Tol 1457-262B 


9.00 


0, 


.792 


0.582 


0, 


,020 


2 


.417 


1, 


,41E+08 


Tol 1457-262C 


8.88 


0, 


.669 


0.359 


0, 


,036 


2 


.099 


7, 


,16E+07 


ESO 566-8 


5.17 


0, 


.505 


0.008 


0, 


,414 


0, 


.693 


3, 


,19E+07 


NGC 5253 A 


9.20 


0, 


.851 


0.756 


0, 


,102 


1, 


,754 


6, 


,82E+07 


NGC 5253 B 


9.38 


0, 


.856 


0.775 


0, 


,086 


1, 


,841 


7, 


,llE+07 


NGC 5253 C 


8.03 


0, 


.773 


0.532 


0, 


,041 


2, 


,046 


2, 


,60E+08 


NGC 5253 D 


7.67 


0, 


.527 


0.048 


0, 


,079 


1, 


,582 


7, 


,72E+07 



" Value derived for the q parameter (in units of cm s ^) obtained using the optimal calibration given bv lKewlev fc Dopital ||2002|') . 
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Table A. 2. Results of the oxygen abundance, in the form 12-|-log(0/H), for objects with a direct estimation of the metallicity, 
considering several empirical calibrations. 



Region Branch Tg EP84 MRS85 M91 ZKH94 POO POl" PTOS*" KD02 KK04 D02 PP04a PP04b PP04c N06a= N06b N06c 
R23 R23 R23,y R23 R23 R23, PR23, PR23, aR23, 9 N2 N2 N2 O3N2 R23 N2 O3N2 



HOG ,il AC 


H 


8 


22±0 


05 


8 


25 


8 


89 


8.67 


8.74 


8.47 


8 


15 


8 


09 


7 


99 


8 


12 


8 


40 


8.34 


8 


29 


8 


30 


8 


05 


8.16 


8 


22 


HCG 31 B 


H 


8 


14+0 


08 


8 


14 


8 


48 


8.29 


8.44 


8.24 


8 


22 


8 


12 


8 


41 


8 


44 


8 


39 


8.33 


8 


28 


8 


29 


8 


07 


8.16 


8 


20 


HCG 31 E 


H 


8 


13+0 


09 


8 


17 


8 


62 


8.14 


8.53 


8.31 


8 


18 


8 


13 


8 


19 


8 


32 


8 


35 


8.30 


8 


26 


8 


25 


8 


07 


8.11 


8 


15 


HCG 31 1 1 


L 


8 


07+0 


06 


8 


02 


8 


30 


8.13 




7.86 


8 


12 


7 


99 


8 


46 


8 


33 


8 


05 


8.07 


8 


09 


8 


03 






7.81 


7 


67 


HCG 31 F2 


L 


8 


03+0 


10 


7 


90 


8 


54 


8.06 


8.48 


7.76 


8 


13 


7 


95 


8 


19 


8 


27 


8 


06 


8.07 


8 


10 


8 


07 


8 


07 


7.83 


7 


79 


HCG 31 G 


H 


8 


15+0 


08 


8 


13 


8 


43 


8.26 


8.40 


8.22 


8 


11 


8 


17 


8 


31 


8 


42 


8 


41 


8.34 


8 


29 


8 


26 


8 


07 


8.17 


8 


16 


Mkn 1199 C 


H 


8 


75+0 


12 


9 


37 


9 


26 


9.00 


9.18 


9.19 


8 


71 


8 


54 


9 


14 


9 


14 


8 


92 


8.74 


8 


90 


8 


81 


9 


18 


8.78 


8 


94 


Haro 15 A 


H 


8 


10+0 


06 


8 


08 






8.14 




7.91 


8 


12 


8 


12 


8 


48 


8 


34 


7 


98 


8.01 


8 


05 


7 


97 






7.74 


7 


38 


Mkn 5 Al 


L 


8 


07+0 


07 


7 


90 


8 


54 


8.04 


8.48 


7.76 


8 


13 


8 


13 


8 


19 


8 


26 


8 


18 


8.17 


8 


16 


8 


12 


8 


07 


7.94 


7 


89 


Mkn 5 A2 


L 


8 


08+0 


07 


7 


95 


8 


43 


8.14 


8.41 


7.81 


7 


92 


8 


17 


8 


18 


8 


33 


8 


16 


8.15 


8 


15 


8 


11 


8 


07 


7.92 


7 


87 


IRAS 08208+2816 


H 


8 


33+0 


08 


8 


15 


8 


50 


8.55 


8.46 


8.25 


8 


42 


8 


35 


8 


35 


8 


25 


8 


47 


8.39 


8 


34 


8 


23 


8 


35 


8.23 


8 


11 


POX 4 


L 


8 


03+0 


04 


8 


15 






8.20 




7.97 


7 


92 


8 


06 


8 


48 


8 


40 


7 


78 


7.86 


7 


91 


7 


87 






7.53 






UM 420 


L 


7 


95+0 


05 


7 


78 


8 


73 


7.98 


8.61 


7.66 


7 


85 


7 


86 


8 


02 


8 


16 


8 


20 


8.39 


8 


33 


8 


28 


7 


57 


8.15 


8 


14 


SBS 0926+606A 


L 


7 


94+0 


08 


7 


88 


8 


57 


7.97 


8.50 


7.75 


7 


77 


7 


80 


8 


17 


8 


20 


7 


97 


8.00 


8 


05 


8 


02 


7 


71 


7.73 


7 


64 


sBb 0948+532 


L 


8 


03+0 


05 


8 


01 


8 


31 


8.06 




7.86 


7 


82 


8 


10 


8 


34 


8 


28 


7 


91 


7.95 


8 


01 


7 


95 






8.01 






SBS 1054+365 


L 


8 


00+0 


07 


8 


05 


8 


21 


8.09 




7.89 


7 


84 


7 


91 


8 


48 


8 


30 


7 


87 


7.93 


7 


98 


7 


93 






7.63 






SBS 1211+540 


L 


7 


65+0 


04 


7 


86 


8 


60 


7.85 


8.52 


7.73 


7 


68 


7 


65 


8 


02 


8 


10 


7 


58 


7.70 


7 


69 


7 


84 


7 


68 


7.31 






SBS 1319+579A 


L 


8 


05+0 


06 


8 


09 






8.13 




7.93 


8 


11 


8 


11 


8 


48 


8 


33 


7 


77 


7.85 


7 


90 


7 


88 






7.52 






SBS 1319+579B 


L 


8 


12+0 


10 


7 


85 


8 


62 


8.01 


8.53 


7.72 


8 


13 


8 


12 


8 


13 


8 


23 


8 


14 


8.14 


8 


14 


8 


11 


8 


07 


7.90 


7 


89 


SBS 1319+579 C 


L 


8 


15+0 


07 


7 


85 


8 


62 


8.02 


8.53 


7.72 


8 


13 


8 


13 


8 


12 


8 


23 


8 


18 


8.17 


8 


16 


8 


13 


8 


06 


7.94 


7 


93 


SBS 1415+437 C 


L 


7 


58+0 


05 


7 


63 


8 


91 


7.72 


8.76 


7.53 


7 


57 


7 


55 


7 


86 


7 


99 


7 


79 


7.86 


7 


92 


7 


99 


7 


39 


7.55 


7 


55 


SBS 1415+437 A 


L 


7 


61+0 


06 


7 


58 


8 


96 


7.64 


8.80 


7.49 


7 


50 


7 


48 


7 


82 


7 


92 


7 


72 


7.81 


7 


86 


7 


97 


7 


34 


7.48 


7 


41 


III Zw 107 


H 


8 


23+0 


09 


8 


17 


8 


62 


8.57 


8.53 


8.31 


8 


40 


8 


35 


8 


13 


8 


24 


8 


39 


8.33 


8 


28 


8 


23 


8 


46 


8.15 


8 


10 


Tol 9 INT 


H 


8 


58+0 


15 


8 


30 


9 


00 


8.76 


8.84 


8.58 


8 


61 


8 


55 


8 


95 


8 


90 


8 


68 


8.56 


8 


54 


8 


42 


8 


77 


8.46 


8 


40 


Tol 9 NOT 


H 


8 


55+0 


16 


8 


29 


8 


97 


8.73 


8.81 


8.55 


8 


56 


8 


50 


8 


94 


8 


88 


8 


72 


8.59 


8 


59 


8 


44 


8 


75 


8.51 


8 


44 


Tol 1457-262A 


L 


8 


05+0 


07 


8 


09 






8.26 




7.92 


8 


11 


8 


20 


8 


58 


8 


42 


8 


04 


8.06 


8 


09 


8 


02 






7.80 


7 


63 


Tol 1457-262B 


L 


7 


88+0 


07 


8 


55 






8.21 




7.87 


7 


91 


8 


21 


8 


38 


8 


29 


7 


88 


7.93 


7 


99 


7 


96 






7.60 






Tol 1457-262C 


L 


8 


06+0 


11 


8 


48 






8.21 




7.88 


8 


00 


8 


24 


8 


48 


8 


37 


8 


07 


8.08 


8 


10 


8 


06 






7.83 


7 


77 


ESO 566-8 


H 


8 


46+0 


11 


8 


27 


8 


92 


8.68 


8.77 


8.50 


8 


44 


8 


38 


8 


92 


8 


84 


8 


84 


8.68 


8 


76 


8 


51 


8 


70 


8.66 


8 


53 


NGC 5253 A 


H 


8 


18+0 


04 


8 


09 


8 


24 


8.13 




8.15 


8 


11 


8 


13 


8 


53 


8 


33 


8 


40 


8.34 


8 


28 


8 


17 






8.16 


8 


00 


NGC 5253 B 


H 


8 


19+0 


04 


8 


09 


8 


21 


8.14 




8.14 


8 


11 


8 


13 


8 


48 


8 


34 


8 


34 


8.29 


8 


25 


8 


14 






8.10 


7 


94 


NGC 5253 C 


L 


8 


28+0 


04 


8 


14 


8 


46 


8.53 


8.42 


8.23 


8 


38 


8 


32 


8 


67 


8 


63 


8 


11 


8.11 


8 


13 


8 


08 


8 


30 


7.87 


7 


80 


NGC 5253 D 


L 


8 


31+0 


07 


8 


15 


8 


52 


8.19 


8.47 


8.26 


8 


23 


8 


17 


8 


32 


8 


37 


8 


31 


8.27 


8 


23 


8 


22 


8 


37 


8.07 


8 


10 



NOTE: The empirical calibrations and the parameters used for each of them arc: EP84: Edmunds & Pagel (1984) that involves the Ri^ 
parameter; MRS85: McCall, Rybski & Shields (1985) using i?23; M91: McGaugh (1991) using K23 and {/; ZKH94: Zaritzky, Kennicutt & 
Huchra (1994) using K23; POO: Pilyugin (2000) using i?23; POl: Pilyugin (2001a,b) using i?23 and P; KD02: Kewley & Dopita (2002) using 
i?23 and q\ KK04: Kobulnicky & Kcwlcy (2004) using R23 & q; D02: Denicolo, Tcrlevich & Terlcvich (2002) using the N2 parameter; PP04; 
Pettini & Pagel (2004), using (a) ^2 with a linear fit, (b) N2 with a cubic fit, (c) the O-iN^ parameter; N06: Nagao et al. (2006) using the 
cubic relations involving the i?23 (a), A'^2 (b) and OiN'i (c) parameters. The value compiled in the column labeled Tc is the oxygen 
abundance derived by the direct method. 

" The value listed for POl is the average value between the high- and the low-metallicity branches for objects with 7.90<12+log(O/H)<8.20. 
^ The value listed for PT05 is the average value between the high- and the low-metallicity branches for objects with 
8.05<12+log(O/H)<8.20. 

° The value listed for N06 is the average value between the high- and the low-metallicity branches for objects with 8.00<12+log(O/H)<8.15. 



30 



